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Abstract 

We study QCD radiation in decay processes involving heavy particles. As 
input, the first-order gluon emission rate is calculated in a number of re- 
actions, and comparisons of the energy flow patterns show a non-negligible 
process dependence. To proceed further, the QCD parton shower language 
offers a convenient approach to include multi-gluon emission effects, and 
to describe exclusive event properties. An existing shower algorithm is ex- 
tended to take into account the pro cess- dependent mass, spin and parity 
effects, as given by the matrix element calculations. This allows an im- 
proved description of multiple gluon emission effects off b and t quarks, and 
also off nonstandard particles like squarks and gluinos. Phenomenological 
applications are presented for bottom production at LEP, Higgs particle 
decay to heavy flavours, top production and decay at linear colliders, and 
some simple supersymmetric processes. 
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1 Introduction 



Heavy particles often tend to be in the focus of particle physics research. The search 
for new physics at the high-energy frontier is an obvious example. Not only are the 
hypothetical new heavy states here of interest, but often their detection relies on decay 
chains that involve heavy flavours. Here the b quark is the prime example. By secondary 
vertices or semileptonic decays it can be tagged. It is a main decay product of top, where 
the b tag is a standard requirement for top identification above background at hadron 
colliders. A light Standard Model Higgs (m h < 140 GeV) predominantly decays to bb, 
and for any Higgs state the bb branching ratio is a key parameter in pinning down its 
nature. These tagging and identification aspects also appear in other scenarios beyond 
the standard model. 

Many of the potentially new particles carry colour, like squarks, gluinos and lepto- 
quarks. Also colourless states often decay to quarks. Therefore QCD radiation is un- 
avoidable, and the large value of the a s coupling implies that such radiation can be quite 
profuse. It is well-known, e.g. from LEP, that fixed-order perturbation theory fails to 
describe QCD radiation effects off light quarks: the rate of a few well separated jets 
can be reproduced, but the internal structure of these jets requires the resummation of 
multiple-gluon effects. The most successful method for achieving such a resummation is to 
apply the parton-shower language, where explicit final states can be generated, with full 
respect for energy-momentum conservation and other constraints. By introducing some 
low fixed cut-off scale Qq ~ 1 GeV, the very soft and nonperturbative regime of QCD can 
be factored off and put in a universal hadronization description, such as the Lund string 
fragmentation model [1]]. With this approach, it is possible to obtain a quite accurate 
description of essentially all hadronic final-state properties of e + e~ annihilation events. 
The physics is more complicated in a hadron collider environment, with additional effects 
e.g. from initial-state QCD radiation and underlying events, so the level of ambition may 
have to be set accordingly. However, this applies to all conceivable descriptions, so again 
we expect a realistic description of jets to be most easily achieved with the help of the 
parton-shower language. 

Heavy quarks radiate less than light ones in the soft-gluon region, so in some respects 
QCD multiple-emission corrections are not as crucial. The b is sufficiently light that 
it still radiates profusely, but top and supersymmetric particles may be heavy enough 
that multiple-gluon emission effects are limited. Furthermore, for particles that are very 
short-lived, the width will provide a cut-off on soft-gluon emission [[J . So for new-particle 
searches, higher-order QCD effects may be less of an issue. However, whenever higher 
precision is required, like for a mass determination of a new state, it may still be necessary 
to model whatever soft-gluon emission effects there are. And, of course, the decays of these 
new states will bring further radiation, e.g. like t — > bW + , where the radiation off the 
t and off the b become intertwined, and where width effects may become important for 
short-lived particles. The aim of the current article is precisely to improve the description 
of gluon radiation off b and c quarks as well as off heavier objects, in order to allow 
higher-precision physics studies. 

The starting point is the calculation of a large number of first-order matrix elements, 
for gluon emission in decay processes, often generalizing on results already found in the lit- 
erature. The processes thereby covered include gluon emission e.g. in 7*/Z°/h°/A — > bb 
(with differences between a vector, axial vector, scalar and pseudoscalar source of the 
b's), t — > bW + , and various supersymmetric reactions. Quite apart from the implemen- 
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tation in the shower framework, these calculations can also be used to assess the degree 
of (non) universality, i.e. the dependence on the colour and spin structure of the processes 
under identical kinematical conditions. In particular, we will show that the 'dead cone' 
picture can be quite misleading for the emission of more energetic gluons. 

These matrix elements are taken as input to an improved version of an existing shower 
algorithm ||, which primarily has been used for radiation off lighter quarks and gluons. 
This shower algorithm tends to slightly overestimate the amount of gluon radiation for 
e + e~ — > 7*/Z° — > qqg when q is massless, so a simple rejection procedure can be used to 
match the hardest emission to the massless matrix elements. The radiation off heavier 
quarks is underestimated in the collinear region, however. Or, put differently, the dead 
cone [HI effects ggerated. In order to allow a corresponding approach for heavy 

flavours, the evolution variable therefore has been modified to bring up the emission rate. 
It therefore becomes possible to use the process-specific matrix elements as rejection 
factors also here, and furthermore to do it in all stages of the shower. 

The one place where precision tests are possible is in studies of b production at LEP1 
f|. We therefore study this topic in detail here, and use it as a test bed for different 
possible variations of the basic shower algorithm. Other studies include Higgs, top and 
super symmetric particle production and decay. Furthermore, we also include some plots 
to illustrate differences between gluon radiation off sources of different spin. 

While not the main thrust of this article, we note that new data on gluon branching 
to heavy flavours, g — > cc and g — > bb, has been presented by the LEP collaborations 
0. We therefore also discuss what implications this might have for the shower algorithm, 
and whether the data could be accommodated by reasonable modifications. 

The plan of the article is as follows. In Section 2, the existing older shower algorithm 
is explained, together with a few intermediate variants thereof. The new approach is 
described in Section 3. A survey of the matrix-element calculations and the consequent 
radiation patterns is given in Section 4. Some applications are then presented in Section 
5. The separate topic of gluon splitting is discussed in Section 6. Finally, a summary and 
outlook is given in Section 7. 



2 Previous Models 

Several shower algorithms have been proposed in the literature. Today the three most 
commonly used ones probably are those found in Pythia/Jetset [[| ||, Herwig |7|, ^ 
and Ariadne ||. The studies in this article will be based on the former one. 

The Pythia final-state shower consists of an evolution in the squared mass m 2 of a 
parton. That is, emissions are ordered in decreasing mass of the radiating parton, and 
the Sudakov form factor ]TD[ is defined as the no-emission rate in the relevant mass range. 



Such a choice is not as sophisticated as the angular one in Herwig or the transverse 
momentum one in Ariadne, but usually the three tend to give similar results for e + e~ 
annihilation events. An exception, where small but significant differences were found, is 



the emission of photons in the shower |TT| . In general, comparisons between the three are 
helpful in estimating a range of theoretical uncertainty, in interpretations of existing data 
or in predictions for the future. 

One of the advantages of the Pythia algorithm is that a mapping between the 
parton-shower and matrix-element variables is rather straightforward to 0(a s ) for mass- 
less quarks, and that already the basic shower populates the full phase space region very 
closely the same way as the matrix element. It is therefore possible to introduce a simple 
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correction to the shower to bring the two into agreement. Also in Ariadne the massless 
matrix-element matching is straightforward. By contrast, the Herwig angular-ordered 
approach does not automatically cover the full qqg phase space, which means that a 
subset of three-jet events are not generated at all. This leads to problems in the descrip- 
tion of LEP data, that were overcome by separately adding the missing class of three-jet 
events ]12| , in addition to matching to the matrix elements in the allowed region. More 



recently, a similar approach has been applied to top decay fl3|] , where again Herwig did 
not populate the full phase space. 

2.1 The massless shower 

In addition to mass, the other main variable in the Pythia shower is z, as used in the 
splitting kernels. It is defined as the energy fraction taken by the first daughter in the CM 
frame of the event. That is, in a branching a — > b + c, Eb = zE a and E c = (l — z)E a . In the 
original choice of z, which is done at the same time as m a is selected, the b and c masses 
are not yet known, and are therefore imagined massless, also in cases where either of them 
is known to have a non-vanishing on-shell mass. A cut-off scale m m \ n = Qo ~ 1 GeV is 
used to constrain the allowed phase space, so that only branchings with m a > m m ; n are 
allowed. For a massive quark, the cut-off is shifted to 



• i Qo , Qo /i \ 

m aMn = \jm? a + — + — . (1) 

The allowed z range, z^ < z < z + , then becomes 

z± = ^ {1 ± f3 a <S>(m a - m aMn )} , (2) 

with P a = \p a \/E a the a velocity and Q(x) the step function. 

At a later stage of the evolution, when nib and m c are being selected, possibly well 
above Qo, the previously found z may be incompatible with these masses. The adopted 
solution is to take into account mass effects by reducing the magnitude of the three- 
momenta pb — —Pc m the rest frame of a. Expressed in terms of four-momenta in an 
arbitrary frame, this is equivalent to 

(°) J- b T,(°) 



p b = (1 - k b )pl + k c p c 

Pc = (1 - k e )p<® + hp® , (3) 

where p^ and p^ are the original massless momenta and pb and p c the modified massive 
ones. The parameters kb and k c are found from the constraints pi = ml and p\ = m^: 

. _ m l ~ X abc ± (ml - m 2 b ) 
b ' c ~ 2ml ' 

with 

X abc = J (ml - ml- m 2 c ) 2 - Amlm 2 c . (5) 
The relation between the preliminary and final energy sharing thus is given by 

z - --(l-k b )z + k c (l-z)- ^— + —z, (6) 
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with z = / E a as above. The transverse momentum p± of b and c with respect to the 
a direction is given by 

but is approximated by ~ z(l — z)m 2 a when used as argument in a s (p 2 ± ) for the shower. 

Angular ordering is not automatic, but is implemented by vetoing emissions that don't 
correspond to decreasing opening angles. The opening angle of a branching a — > b + c is 
calculated approximately as 



'zfl - zl E, 



(8) 



The procedure thus is the following, Fig. |I]a. In the 7*/Z° decay, the two original 
partons 1 and 2 are produced, back-to-back in the rest frame of the pair. In a first step, 
they are evolved downwards from a maximal mass equal to the CM energy, with the 
restriction that the two masses together should be below this CM energy. When the 
two branchings are found, they define mi and m 2 and the z values of 1 — > 3 + 4 and 
2^5 + 6. These branchings obviously have smaller opening angles than the 180° one 
between 1 and 2, so no angular-ordering constraints appear here. A matching procedure 
to the matrix element is used to correct the branchings, however, as will be described 
below. In subsequent steps, a pair of partons like 3 and 4 are evolved in parallel, from 
maximum masses given by the smaller of the mother (1) mass and the respective daughter 
(3 or 4) energy. Here angular ordering restricts the region of allowed z values in their 
branchings, but there are no matrix-element corrections. Once and are fixed, the 
kinematics of the 1 — > 3 + 4 branching needs to be modified according to eq. This is 
the reason why the evolution is always done for a pair of partons (whereof not both need 
branch further, however), and why the final kinematics of a branching is postponed to a 
later stage than the choice of z value. 

Several other aspects of the shower could be discussed, such as the choice of non- 
isotropic azimuthal angles to improve the coherence description and include gluon spin 
effects, the possibility also to emit photons, or the option to force some branchings in 
order to match higher-order matrix elements fl4] . These are of less interest here and thus 
not covered. 
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2.2 The massless matrix element correction 



Let us now compare the parton-shower (PS) population of three-jet phase space with the 
matrix-element (ME) one, for the case of e + e~ annihilation to massless quarks. With the 
conventional ME numbering q(l) q(2) g(3), with Xj = 2Ej/EcM, the matrix element is of 



the form [15| 



J_ da M E = «s c ^ x\ + x\ 
a$ dx\dx 2 2tc (1 — Xi)(l — X2) ' 

with the colour factor Cp = 4/3. We have normalized to the lowest-order cross section 
cr , so that the expression can be interpreted as a probability distribution. For future 
reference, we will use the notation 

1 da a s N(x u x 2 ,r) 

^ftz rrz r j l-LUJ 



do dx\ dx2 2tx (1 — x±)(l — X2) 

where r = mjE C u = m^/ E CM . Thus N ME (xi,x 2 , 0) = x\ + x\. 

There are two shower histories that could give such a three-jet event. One is 
7*/Z°(0) — > q(i)q(2) — > q(l)q(2)g(3), i.e. with an intermediate (2) quark branching 
q(i) — > q(l)g(3), illustrated in Fig. [l|b. This gives 

Q 2 = m* = (p -p 2 ) 2 = (l-x 2 )E 2 CM , (11) 
PoP ± E x xx xi 



PoPi Ei x x + x 3 2 - x 2 ' 
dQ 2 dz dx 2 dx\ 
Q 2 1 - z 1 - x 2 x 3 

The parton-shower probability for such a branching is 

1 — xi dxi dx 2 



(12) 
(13) 



a s 1 + z 2 dQ 2 _ol s 
2tt 1 — z Q z 2tt 



1 



%2 



x 3 (l—xx)(l—x 2 ) 



(14) 



In the second shower history, the roles of q and q are interchanged, i.e. x\ <-> x 2 . This 
is the same set of Feynman graphs as in the matrix-element description, except that the 
shower does not include any interference between the two diagrams. The two shower 
expressions can therefore be added to give the overall shower population of the three-jet 
phase space, of the form in eq. fllPP but with 



N PS (^x 2 ,0) = ^ 1+ +— 1+ b^M , (15) 




In spite of iVps being lengthier than A^me, it turns out that the two almost exactly 
agree over the whole phase space, but with the shower rate somewhat above, see e.g. 
Fig. 0. It is therefore straightforward and efficient to use the ratio 



dq-ME _ N ME (x 1 ,x 2 , 0) 
d(T PS N PS {xi,x 2 ,0) 



#me/ps(zi, X 2 , 0) = — = Ar z nN (16) 



as an acceptance factor inside the shower evolution, in order to correct the first emission 
of the quark and antiquark to give a sum in agreement with the matrix element. 

Clearly, the shower will contain further branchings that modify the simple result, e.g. 
by the emission both from the q and the q, but these effects are formally of 0(a 2 ) and thus 
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Figure 2: The gluon emission rate as a function of emission angle #13 = # qg , for a 
10 GeV gluon energy at Eqm = 91 GeV, and with rrib = 4.8 GeV. All curves are 
normalized to the massless matrix-element expression, eq. (||), here thus represented 
by the small-dotted line at unity. Dashed: the massless shower before correction, 
Nps(xi, x 2 , 0)/iV ME (xi, x 2 , 0) = l/R ME / PS (xi, x 2 , 0). Full: the rate from massive ma- 
trix elements, N ME (xi, x 2 , r)/N ME (xi, x 2 , 0). Dash-dotted: the rate from massive parton 
shower, iVp S (a;i, x 2 , r)/N ME (xi, x 2 , 0). Large dots: the new rate from massive parton 
showers, A^ s (^i, x 2 , r)/N ME (x 1 ,x 2 , 0). 

beyond the accuracy we strive to match. One should also note that the shower modifies 
the distribution in three-jet phase space by the appearance of Sudakov form factors, and 
by using a running a s (p 2 ± ) rather than a fixed one. In both these respects, however, the 
shower should be an improvement over the fixed-order result. 

2.3 The massive matrix element correction 

The prescription of correcting the first branchings by the factor in eq. (|T6| ) was the original 
one, used up until Jetset 7.3, for massless and massive quarks alike, where the Xi = 
2Ei/EcM variables were defined with masses included. In version 7.4 an intermediate 
'improvement' was introduced, in that the massive matrix element expression was used, 
as given for a vector source like the 7* : 

N M b(xi,x 2 , r) = x\ + x 2 2 - Ar 2 x 3 - 8r 4 - (2r 2 + 4r 4 ) f ^— — + jj— — ) . (17) 

\1 - X\ 1 — x 2 J 

The shower algorithm itself was not changed, nor the assumed shower weight, i.e. an ac- 
ceptance factor N ME (xi, x 2 , r)/Np S (xi, x 2 , 0) was applied for the first branching on either 
side. (The older behaviour remained as an option.) 

The mass suppression in the matrix element is illustrated in Fig. 0. We remind that, 
in the soft-gluon limit, a spin-independent (and thus universal) eikonal expression holds 
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d(J qqg / ,\ ( Pi P2 V d 3 p 3 



oc -1 



P1P3 P2P3J E 3 

2 



/ 2piPa ml 2 } E dE d cos e 

\(PlP3)(P2P3) (PlPs) 2 (P2P3) 2 J 3 3 

In the limit of small angles #13 this gives a mass suppression factor 



da{x 3 ,6 13 ,r) ( 9 



da(x 3 ,9 13 ,0) \9 2 3 + 4r 2 



; 13 



(19) 



i.e. the characteristic dead cone of opening angle approximately 2r = m q /E q . Note that 
the mapping between (xi,X2) and (x 3 , 6*13) depends on r, so eq. ( |19|) is not quite the same 
as iV ME (xi, x 2 , r)/N ME (xi, x 2 , 0). 



2.4 The massive phase space correction in the shower 



More recently ||, the issue of masses in the shower was further studied, since the expres- 
sion for Np S (xi,x 2 ,0) had not been touched when N ME wa s improved. 
In the derivation of Nps(xi, x 2 , r), one can start from the ansatz 



X 2 



x 3 



2 2 

mf — mz 



1 + 



E 2 

2 

mf 



E 2 



mf 



mz 



E 2 



{{l-k x )z + h{l-z)) , 
((1 - k 3 ){l - z) + hz) . 



(20) 



The quark mass enters both in the energy sharing between the intermediate quark i and 
the antiquark 2, and in the correction procedure of eq. (|3[) for the splitting of energy in 
the branching q(i) q(l)g(3). The constraints p\ = m 2 and p\ = give k\ — and 
k 3 = m 2 Jm 2 . One then obtains 

Q 2 = m 2 = (: 



1 



x 2 + r )E CM , 
2 x 3 



2-x 2 
dQ 2 dz 



X\ — r 



l — x 2 
dx 2 



dxi 



x 2 + r 2 x 3 



(21) 
(22) 

(23) 



This gives the answer 
N PS (x 1 ,x 2 ,r) 



X\ 



x 2 



x 3 

+ {xi <- 



1 — x 2 + r 2 
x 2 } , 



1 + 



(2 - x 2 ) 



X\ — r 



x 3 



x 2 



(24) 



where the second term comes from the graph where the antiquark radiates. 

The mass effects go in the 'right' direction, Nps(x\, x 2 ,r) < Np$(xi,x 2 , 0), but actually 
so much so that Nps(xi, x 2 , r) < NuE(xi,X2,r) in major regions of phase space. This is 
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illustrated in Fig. 0. Very crudely, one could say that the massive shower exaggerates the 
angle of the dead cone by about a factor of two, in this rather typical example. There is 
no dead cone as such built in, however, but a rather more coincidental mass suppression 
mainly generated by the factor (1 — x 2 ) /(I — ^2 + r 2 ). 

Thus the amount of gluon emission off massive quarks is underestimated already in 
the original prescription, where masses entered in the kinematics but not in the ME/PS 
correction factor. When the intermediate 'correction' ratio A^me(^I) ^2, r)/Nps(xi, X2, 0) 
is applied, the net result is a distribution even more off from the correct one, by a factor 
Np$(xi, X2, r)/Nps(xi, X2, 0). Thus it would have been better not to introduce the mass 
correction in JETSET 7.4. 

Based on the results above, one can now instead use the correct ME/PS factor 
Nme(xi, X2, r)/Nps(xi, X2, r). A technical problem is that this ratio can exceed unity, 
in the example of Fig. ||] by up to almost a factor of two. This could be solved e.g. by 
enhancing the raw rate of emissions by this factor. However, another trick was applied, 
based on the facts that the shift of eqs. (|) and (|20"D implies that smaller-energy gluons 
would be allowed for a massive quark than a massless one, and additionally that the acces- 
sible z range is overestimated in the original ansatz. Therefore, without any (noticeable) 
loss of phase space, z can be rescaled to a z' according to 



The ME/PS correction factor then has to be compensated by k, and thereby comes below 
unity almost everywhere — the remaining weighting errors are too small to be relevant. 

This procedure, default since Pythia version 6.130, improves the shower description of 
mass effects in the amount of three-jet events ||. Mass effects are only included correctly 
for the first branching of the q and q in the shower, however. Subsequent emissions involve 
no correction procedure. Instead the dead cone effect is exaggerated, similarly to what 
was shown in Fig. [| Furthermore, even for the first branchings, only the possibility of a 
vector source decaying to two identical-mass quarks is included, while the Z° actually is 
a mixture of vector and axial vector, and e.g. the W ± would decay to two unequal-mass 
quarks. We will therefore next try to develop a more powerful and general approach. 



One of the advantages of the Monte Carlo approach is that, so long as some upper estimate 
can be found that allows simple generation, rejection down to a more complex expression 
is straightforward. In particular, for an evolution in some variable Q with a (Sudakov) 
form factor built up from the no-emission probability, the 'veto algorithm' || can be used, 
in which the initial overestimate of the emission rate is compensated by the possibility of 
a rejection, with a continued evolution downwards in Q 2 from the rejected value onwards. 
As we have seen, the current set of Q 2 and z variables is not so convenient, since the 
emission rate off massive quarks is an underestimate of the correct rate, and therefore the 
standard procedure does not work except after some extra tricks. We will therefore pick 
another set of variables, preferably such that they reduce to the old ones in the massless 
limit. Several approaches could have been taken, but here we have chosen the minimal 
one of retaining the z definition and thereby the existing kinematics machinery. As we will 
show, a modification for the branching a — > be from Q 2 = m 2 a to Q 2 = m 2 — (m 2 ) on -shcii 
is enough to do the job. 




with k 




< 1 . 



(25) 



3 A New Approach 
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3.1 The choice of shower variables 



Again consider e + e~ - 
assume Q 2 = m| — 

to 1/Q 2 being the propagator of the off-shell parton 1. Then eq. ( p2|) is unchanged while 



7*/Z° — > q(l) q(2) g(3), and the formulae of eq. (|20|), but now 
• l - m\ for the q(i) — >■ q(l)g(3) branching. That corresponds 



Q 2 



m 2 ■ 
dQ^_ 

g 2 i 



i 

dz 



dx 2 dx\ 



x 2 x 3 



(26) 
(27) 



So, with this simple trick, we have recovered the Jacobian of the massless case, eq. (|T3| 
This means that there is little mass suppression left in the shower evolution proper: 



Np S (x 1 ,x 2 ,r) 



x 3 



1 + 



(2-x 2 ) 



X\ 



x 3 



1 - x 2 



+ {xi <-> x 2 } , 



(2f 



as can be seen in Fig. ||. In particular, note that N PS (xi, x 2 , r) reduces to Np$(xi, x 2 , 0), 
eq. (n), in the soft limit x 3 —>■ 0. Therefore a matrix-element correction now is even more 
required, but also simpler to implement, to bring down the rate to a reasonable level. 

So far, we have assumed decay to two equal-mass quarks, like in Z° decay. For many 
of the examples to be discussed, like W ± or t decay, the two decay products have unequal 
masses. It is also convenient to allow slightly unequal masses for particles with non- 
negligible widths, e.g. in e + e~ — > tt. We will therefore generalize to the case with 
ri = m l /E C yi = mjE CM ^ r 2 = m 2 /E CM = m^/ 'Ecu- The range of kinematically 



allowed Xj values is then 2r\ < x% < 1+r 2 - 
with the joint condition that 



r 2) 2r 2 < x 2 < l+r^—rf, < x 3 < l — {r\+r 2 ) 



(2 - 2x x - 2x 2 + x x x 2 + 2rf + 2r 2 ) 2 < (xj - Ar{)(x 2 2 - Ar 2 ) . 



The ansatz of eq. (20) is modified by m q — > m 2 while k 3 — > m\/m 2 . Furthermore 



Q 2 



rri; 



(P0-P2) 2 - ml = (1 



x 2 )E t 



2 

CM 5 



1 



2-x 2 
dQ 2 dz 



x 1 



M+r 2 
dxo 



L ~X 2 

dxi 

Q 2 1 — z 1 + r\ — r\ — x 2 x 3 



(29) 

(30) 
(31) 

(32) 



When a colourless particle decays to two colour triplets, the two possible one-gluon- 
emission shower histories then add up to give 



do 



PS 



cr dx\ dx 2 



2« F \x 3 (l + r 2 -r 2 



x 2 



+ 



l + z 2 



x 3 (l + vf - r% - x t ) 



(33) 



Here z\ is given by eq. ( |5T| ) and z 2 is obtained by exchanging 1 <-> 2. Since the numerators 
1 < 1 + z 2 < 2, their exact form is not a main concern for the qualitative discussions. 
As it turns out, the Monte Carlo procedure is simplified if the shower is generated with 
a numerator 2 instead of 1 + z 2 , so we will make this replacement in the following, with 
the matrix element correction procedure compensating this overestimate. 
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3.2 The matrix element correction 



Next, one should try to relate this to the structure of the matrix elements under similar 
conditions. We may expect graphs where either of the partons can radiate a gluon and 
therefore can give a propagator 1/Q\ — l/(m 2 3 — mf), cf. eq. (p0|), or \jQ\ = l/im^— m 2 ,), 
depending on which side radiates. (Diagrams with a four-boson vertex, such as 7*/Z* — > 
ttg, are not singular and therefore do not affect the discussion.) After summing and 
squaring, the cross section then should have the form 

1 do-ME a s f E^ M A(x 1 ,x 2 ) E^ M B(x 1 ,x 2 ) E^ M C(x 1 ,x 2 ) 

s — — o oT^ — I ; — o ^tt; — r 



ctq dxi dx 2 2n \ (m 2 3 — m 2 ) 2 (m 2 ^ — m 2 .) 2 ( m i3 — m i)( m 2 



2:i 



^ r \ A(xi,x 2 ) B(x 1 ,x 2J 
2tt F \(l + r 2 -r 2 -x 2 ) 2 (l+r 2 -r 2 -x^ 2 
C(xi,x 2 



1 + r| — rf — x 2 )(l + r{ — r\ — x\) 



The individual functions A, B and C depend on the gauge choice, but the total cross 
section of course is gauge-independent. It also has some general features, which can be 
seen in the soft-gluon limit expression, eq. (Ti~8|). Here the interference term is giving 
a positive expression and the seemingly quadratically-divergent terms are negative and 
drive the cross section to zero in the collinear limit, the dead cone. In order to have an 
upper estimate of the cross section, for Monte Carlo applications, it is thus enough to 
have the singularity structure of the interference term modelled, and no need to worry 
about quadratic divergences. 

As before, we could now compare the total shower and total matrix element rates, and 
define a corrective factor between the two. To simplify some of the continued studies, 
however, we have adopted an alternative approach. Instead of adding the two shower 
histories to compare with the full matrix element, one can split the matrix element in two, 
such that each part can be compared with only one shower history. Such a subdivision of 
course is arbitrary, but should still be sensible. A gluon emitted close to the pi direction 
should predominantly be emitted by parton 1, and vice versa. A suitable such subdivision 
is in the proportions given by the propagators, 1/Qf '■ ^-/Q 2 , which also is the proportions 
between the two shower histories (in the 1 + z 2 — > 2 approximation). For massless partons, 
and going to the soft-gluon limit, this means a probability (1 + cos 6 l i 3 )/2 = (1 — cos 6 2 ^)/2 
for emission off parton 1 and the rest for emission off parton 2. 

In the description of the gluon emission rate off parton 1, a matrix element fraction 

TJ/ / % Ql 1 do-ME l + rf-rj-X! 1 do- M E , oc , 
IW^H^^-g^- (35) 

should be compared with the first half of the total parton shower expression in eq. (pB]), 

a 2 
W/ps,^, x 2 ) = - CF x3{1 + r2 _ r2 _ x2) • (36) 

The ME/PS correction factor then becomes 

W M E,l(^l,^2) 



Rx(xi,x 2 ) 



W PS ,i(xi,x 2 ) 

(1 + r\ - r\ - x 2 )(l +rf - r\ - x\) ( a s \ - 1 1 da; 



?C F ) -f^, (37) 
2.71 / a dxi dx 2 
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which reduces to Ri(xi, x 2 ) = N M - E (xi,x 2 ,r)/2 for r x = r 2 = r. 

The intention is that Ri(x\,x 2 ) should be finite and well-behaved over all of phase 
space, with one factor of each divergence now multiplied on to the matrix element ex- 
pression. We will illustrate the behaviour later on, but the key observation is that, for all 
the matrix elements we will study, searches over the full mass parameter plane and phase 
space have failed to find any point where the ratio is above unity. R\ therefore serves well 
as simple rejection factor. (Without the replacement 1 + z 2 — > 2 in the shower description 
one would find Ri(xx,x 2 ) > 1, and some extra precautions would be necessary.) The 
correction factor R 2 for the other shower history looks the same. Put another way, the 
relative probability for a gluon to be emitted by parton 1 or by parton 2 is unchanged by 
the matrix element correction procedure. 

The above procedure works not only for quarks, but also e.g. for squarks. The 
numerator of the splitting kernel would have been 2z rather than 1 + z 2 |Tj|, but both 
are equally well approximated from above by 2. For a gluino, the colour charge is Nc 
rather than Cf, so the assumed shower emission rate has to be scaled up by a factor 
Nc/Cf = 9/4 (also for the recoiling colour triplet parton, since the separation of radiation 
is not perfect), with no other change required. 

We will also encounter processes where the decaying particle carries colour, like in t — > 
bW + . Then gluon emission off this particle has to be considered, i.e. graphs like t(0) — > 
t(z) g(3) — > b(l) W + (2) g(3), proceeding via an intermediate off-shell top. Neglecting the 
top width, this introduces a new kind of inverse propagators 

IQoI = m o - m i = m i23 - m 12 = 2 (Pl + P?)P3 = x 3 E^M . (38) 

The leading term of the matrix element, proportional to 1/(QqQ1), will therefore have 
the same singularity structure as the shower rate in eq. fl36|). In a shower description, 
where only the b is allowed to radiate, there is thus no problem in principle of letting that 
radiation account for the full emission pattern of the matrix elements, i.e. 

W M e(x 1 ,x 2 ) 



Ri(x 1 ,x 2 ) 



Wps,i(zi,£2) 

x 3 (l + r\ - r\ - x 2 ) (a s X- 1 1 dcr 



(— Cf) — -7—7^- ■ (39) 
\2n ) (in dxi axo 



2 \2n J cr dxi dx 2 

Since also here the shower rate turns out to exceed the matrix element one, a simple 
rejection approach should work well. 

3.3 Subsequent branchings 

So far, we have considered matrix elements as providing the probability of exclusive 
three-jet events. An alternative interpretation, however, would be in terms of an inclusive 
density of gluon emissions, with the possibility of several such per event. This interpre- 
tation works well in the soft-gluon limit, while the emission of a hard gluon reduces the 
phase space for subsequent emissions and thus ruins the picture of independent emissions. 
Furthermore, the possibility for gluons to branch in their turn leads to the need to in- 
clude coherence effects [T7) that constrain allowed emissions. Nevertheless, the matrix 



elements can be used to extract important information, not only on the emission rate of 
hard gluons, but also on that of soft and collinear ones. 

Historically, different approaches have therefore been taken. In the Jetset/Pythia 
procedure used until now, only the first branching is corrected by matrix-element infor- 
mation. Subsequent emissions involve no corrections, but only the process-independent 
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splitting kernels, and therefore give the wrong emission rate off heavy quarks that has 
already been noted. In the Herwig routine, a correction is performed for every emission 
that is the hardest so far JT2| - whereas Jetset emissions tend to be ordered in hard- 



ness, this is less likely in Herwig, so the distinction is then relevant. Ariadne, finally, 
imposes a correction at all steps of the cascade. 

Since our new algorithm does not have the correct dead-cone behaviour built in from 
the onset, it is clear that some correction procedure will be required. Rather than im- 
posing the process-independent collinear behaviour, we have chosen to base ourselves on 
the matrix elements also in this region. This offers a smooth interpolation between the 
hard process-specific and the soft collinear universal behaviours, and sidesteps the issue 



of which emissions are to be considered the hardest so far [12| . Thus the shower rate will 
be corrected to the matrix element one at every step of the shower off the original par- 
tons. The gluon cascading g — > gg is unaffected, since there are no gluon-mass effects or 
matrix elements to be considered here. Light quarks (in 7*/Z° decay) are also essentially 
unaffected, since there the shower matches so well with the matrix elements anyway. 

The kinematics of the cascade changes in each emission, as the energy of the radiating 
parton is reduced by the previous emissions. The mapping of an emission on to the matrix 
element variables is thereby not unique. However, from the dead cone formula, eq. fllTf ), 
we see that the emission angle and the mass-to-energy ratio of the emitting quark should 
be represented faithfully in the choice of matching matrix element variables. For the 
branchings subsequent to the first one, say 3 —>■ 7 + 8 in Fig. |l|a, the mother energy E 3 is 
fixed. The choice of a Q 2 = m\ — m 2 - and a z of the branching maps onto four-momenta 
•p-j and ps as described in eqs. (0)^(0) > with p 2 = > and pi = 0. In order to make 
contact with the matrix element variables, now construct a hypothetical recoiling parton 
2', such that p' 2 = — P3 and p' 2 2 = m 2 ^. The CM energy of this reduced system is then 
given by 



E' CM = E 3 + E' 2 = E 3 + ^/p 2 + m* = E 3 + ^ E 2 - Q 2 - m\ + w 2 . (40) 

Now matrix element and parton shower weights can be evaluated for x[ = 2E- 7 /E' CM , 
x 2 = 2E 2 / E' CM , r[ = m q / E' CM and r' 2 = m^/E' CM . Note that, in the limit E± — ► 0, also 
E' CM — > Eqm, and the ordinary matrix elements are recovered. 

The above procedure is not unique. As an alternative, one could retain the original 
Eqm, but construct an off-shell p 2 to carry all the energy and momentum of the system 
except for p 3 = p 7 +p&, i.e. p 2 = p 2 +p 4 in the case of Fig. |TJa. This gives equivalent results 
so long as the gluon energy already emitted is not too large, and else gives a somewhat 
lower rate of wide-angle emission, reflecting that the parton 2' then is assigned such a 
large mass that it radiates less. The overall picture therefore is less appealing, even if 
results in practical applications are almost equivalent. 

In this article, emphasis is put on the gluon emission off the primary quarks (or other 
primary particles). The subsequent branchings g — ► gg have not been affected by the 
mass considerations, so are not discussed here. The rate of g — > qq branchings is another 
topic of some interest, given that LEP results do not quite agree with predictions ||. 
This issue is further discussed in Sect. 6. 

However, given that such a secondary cc or bb pair has been produced, the possibility 
of further gluon emission off this pair should be considered, even if most pairs are produced 
at such a low mass that the phase space left for further radiation is limited. In order to 
provide a sensible behaviour in the collinear region, again matrix-element input is applied, 
calculated for the decay of a colour octet source. (Process 66 of Table [l|, so with the wrong 
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spin of the source, but correct for the radiating parton, which is the main point.) To first 
approximation, this means that radiation occurs independently off the q and q; see the 
discussion on radiation patterns below. The kinematics for the matrix element corrections 
is set up about as described above, i.e. by mapping onto a reduced system at rest with 
preserved energy for the radiating parton. A simplification is that here the recoiling 
parton always has the same mass as the radiating one. 



3.4 Additional issues 

In the older shower algorithm, the decay products of a branching a — > be were assumed 
massless until assigned a mass by some subsequent step. This meant that the p\ of a 
branching, eq. ([?]) with (3 a = 1 as a first simplification, was simplified further to p\ ~ z(l — 
z)m 2 a when used as argument in a s (jr^). In a branching q — > qg, where the q is assumed to 
have a non- vanishing rest mass, one would instead obtain p\ « z(l — z)m 2 (l — ml/m 2 ) 2 . 
(This follows from eq. ([F]) or alternatively from the rescaling of X3 in eq. (pop.) A smaller 
p 2 ± estimate implies a larger a s and emission rate for a given kinematical configuration, 
but also that more phase space is cut out by the requirement p\ > Q§/4, where Qo is the 
soft cut-off scale of the cascade. The net result of such a change is therefore not obvious, 
and we will study it later on. 

Also the calculation of the approximate opening angle of a branching a — > be, eq. (§), 
would be affected by the same considerations. The massive parton energy is increased 
at the expense of E c , as given by eq. (|20|) . Since the common p± is also decreased, the ratio 
P_l/E c is preserved, while p±/E b is decreased. The net result is a decrease of the opening 
angle by a factor (1 + (ml/m 2 )(l — z)/z)~ 1 . As another option, we will then consider 
the consequences of such a decrease in the decay opening angle, without any change of 
the production angle, as a way of minimally relaxing the angular ordering condition off 
massive quarks. 

While our choice of Q 2 ~ m 2 variable has significant advantages for the matching 
to matrix-element expressions, it does not offer as neat an implementation of coherence 
effects as the angular variable of Herwig or the transverse momentum one of Ariadne. 
Without any further constraints, the amount of radiation is overestimated. Therefore, by 
default, angular ordering of emissions is imposed as a further constraint. This, on the 
other hand, tends to restrict emissions somewhat too much. The bulk of these ambiguities 
affect rather soft gluons, which do not give rise to separate jets. For the precision studies in 



Section f>.l\ however, also small effects could be of interest. We have therefore introduced 
a new "intermediate" coherence option, as compared to the minimal modification above, 
wherein no angular constraint is imposed on emissions off the primary qq pair. These 
emissions thus are ordered only in mass. Angular ordering is still imposed in the cascades 
initiated by the gluons emitted off the primary quarks. In particular each such cascade 
is restricted to a cone given by the emission angle of the initiating gluon. In this option, 
generic event properties are only slightly changed compared with the default procedure. 

So far, we have mainly considered configurations where a colour singlet decays to 
stable particles, e.g. 7*/Z° — > bb. Another class of events involve sequential decays of 
coloured objects. The obvious example would be 7*/Z* — > tt — > bW + bW~. Even leaving 
aside the continued fate of the Ws, e.g. assuming they decay leptonically, the event now 
contains four colour charges that may radiate. In the limit T t — > 0, the radiation in the 
top production stage 7*/Z* — > tt decouples completely from that in the top decays. For 
a finite T t , gluons with energies below or around this scale can receive contributions from 
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many colour sources, however, leading to complex radiation patterns ||. Close to the tt 
threshold, the main sources are the respective top decays, described by dipoles tb and 

tb, and radiation off the bb dipole created after the t and t have decayed. Gluons with 
energies above (below) T t predominantly feel the former (latter) dipoles. The radiation 
pattern can be written as 

e 2 /~ -\ r 2 

tb + tb)+-^— bb. (41) 



Thus the reduced radiation induced by the top finite lifetime is compensated by the 
radiation from new sources that would not have been present for a long-lived top. The 

bb dipole introduces a dependence on the opening angle between the b and the b that 
is not there in the separate top decays, so the compensation is not complete. On the 



perturbative level, this gives a dipole effect [[19] that could be observed at low momenta. 



However, even in the approximation of allowing radiation only within each top quark 
decay separately, the string fragmentation picture |l] would imply that a nonperturbative 
colour string should be stretched between the b and b (or their respective cascades), and 



this would introduce a string effect |20| of almost equal character and magnitude pl| . 
Thus only very careful studies, e.g. for high-precision measurements of the top mass, 
would be sensitive to the detailed nature of the soft-gluon emission source. The critical 
transition is the one between a top long-lived enough to produce top hadrons and one 
decaying too rapidly for that, where the hadronic final state does change character. 

The shower algorithm contains options that allows it to be run, either with soft or 
with hard emission dampened according to the respective factors in eq. ([H|). However, 
currently the Pythia program contains no machinery to detect when a description of this 
kind is required, nor a prescription how to combine all possible sources of radiation. This 
is obviously an interesting task for the future, but one that will be required primarily for 
particles with a width significantly larger than that of the top. 

The normal generation sequence therefore contains a set of separated showers. For 
instance, in the top example above, the 7*/Z* — > tt induces a first cascade, whereby 
the t and t shower down to the mass shell. Thereafter follows the respective top decay, 
t — > bW + and t — > bW~, and the separate radiation in those decays. At some yet later 
stage, the two W's may decay to quarks that again radiate. Anytime a colour singlet is 
exchanged, like the W's above, there is a clean separation into disjoint QCD subsystems, 
while the exchange of a coloured state like the top will hook up separately showering 
systems to the strings that later will produce the observable hadrons. Interconnection 
effects could complicate this picture, but so far the evidence is that any such effects 
would be small. 

In hadronic collisions, more complex processes would occur, and also initial-state QCD 
radiation has to be considered as a potential source of further interference effects. In this 
article we will not address these additional complications, but defer that for some future 
study. Currently such interference is almost completely neglected, except for some angular 
restrictions [^3] . 



Even if T t is neglected in the showering, for the sake of providing a unique separation of 
radiation before and after the top propagator, this does not mean that all top quarks have 
to have the same mass. Instead, resonance masses are chosen according to the relevant 
Breit-Wigners, convoluted with the respective cross section formulae. In particular, this 
means that the t and t masses of an event would be unequal. The shower algorithm 



14 



described above then operates on these event-by-event masses, without any reference to 
their nominal on-shell equivalents. 

The hadronization of a partonic configuration, obtained by the chain of decays and 
showers outlined above, is described by the Lund string model |IJ. All coloured partons 
belong to strings, stretched from a quark endpoint via a number of intermediate gluons to 
an antiquark one. (Also diquark endpoints and closed gluon loops can be considered, but 
are of no relevance in this article.) Normally each such string would have a reasonably 
large invariant mass, enough to produce several hadrons. However, occasionally a string 
could come to have a small mass, e.g. by the splitting of a string in two by shower 
branchings g — > qq. Then a special treatment may be required for one- or two-hadron 
decays. While the normal string algorithm has remained essentially unchanged over a 
number of years, this low-mass 'cluster' treatment has recently been improved |24|j . The 
direct consequences for the topics studied here are minimal, except that quark mass 
values have now been optimized, especially to describe charm asymmetries in fixed-target 
experiments. The current default values thus are m u = = 0.33 GeV, m s = 0.5 GeV, 
m c = 1.5 GeV and m b = 4.8 GeV. 



4 Matrix elements 

As input and starting point for our shower studies, we will need the matrix elements 
for the processes of interest. In this article this essentially means the two-body decay 
of a particle, with associated gluon radiation. Some of the formulae are available in the 
literature, but most are not, or at least not easily found. We have therefore calculated a 
number of processes. This also gives us a chance to test the degree of universality of the 
radiation patterns in channels of different colour and spin structure, but with the same 
masses. These results are interesting in their own right. 



4.1 Calculations 



A number of matrix elements have been calculated, using Comphep [[25[ for the ac- 
tual calculation, including an extension package for Supersymmetric processes [J26|, and 



Mathematica [[27]] for subsequent simplification of the expressions. The list of lowest- 
order (LO) processes is given in Table [I]. The LO expression gives the two-body decay 
rate of a particle, i.e. a — > be, and the matching first-order (FO) one the same decay with 
an additional gluon in the final state, a — > beg. 

While the matrix element calculations in this section have been performed from 
scratch, some checks are based on results in the literature. These include: V — > qq 
H° 



for m q = rriq 



ie 



qq for m c 



rria 



28], V — > qq for 



m~ 



181, and t -> bW+ 



for mb = 13 



No doubt, many more are available, without our knowledge. 
The process selection is based on the particle content of the Minimal Supersymmetric 
Standard Model (MSSM), i.e. includes squarks q, gluinos g, neutralinos and charginos x, 
and Higgs states h°, H°, A and H 1 * 1 . The idea, however, is that these situations could 
represent also a number of other non-standard particles. For instance, the decay of a spin 
leptoquark Lq — > q£ is closely similar to q — > q%. 

All calculations have been performed in the zero-width limit of the decaying state 
and the decay products, in order to allow a gauge invariant separation of radiation in 
the production and decay stages. (As explained above, the Monte Carlo simulation of 
processes does include mass selection according to the appropriate Breit-Wigners, so we 
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colour 

3 + 3 
3 + 3 
3 + 1 
3 + 3 
3 + 1 
3 + 3 
3 + 1 
3 + 3 
3 + 1 
3 + 3 
3 + 1 
3 + 1 
3 + 3 
3 + 8 
3 + 8 
8 + 8 



spin 



1 + 1 

1 + 1 

- + 1 

2^2 

| + o 

+ 



0^0 + 1 
0^0 + 
0^0 + 



i + 

i + i 

+ i 



75 

1,75, 1+75 
1,75,1+75 

1,75, 1 + 75 
1,75, 1 + 75 

1 

1 

1 

1 

1,75, 1+75 
1,75,1+75 
1,75,1 ±75 
1,75, 1 ±75 
1,75, 1 ±75 
1,75, 1 ±75 



example 
(eikonal) 
Z° — > qq 
t ^bW+ 
H° -> qq 
t^bH+ 
Z ^qq" 
q -> q'W+ 
H° 



q 



x 
q 
t 

g 
q 
t 



-> qq 
q'H+ 
-> qq 

■» qx 
+ tx 

-> qq 
-> qg 
-tg 



ciKonai 



codes 
6-9 
11 - 14 
16 - 19 
21 - 24 
26 - 29 
31 - 34 
36 - 39 
41-44 
46 - 49 
51 - 54 
56 - 59 
61 - 64 
66 - 69 
71 - 74 
76 - 79 
81 - 84 



Table 1: The processes that have been calculated, also with one extra gluon in the final 
state. Colour is given with 1 for singlet, 3 for triplet and 8 for octet. See the text for an 
explanation of the 75 column and further comments. 



here only comment on the width dependence of the additional gluon radiation, i.e. in 
the ratio of first to leading order cross sections.) Such a separation occurs naturally 
for exchanged colourless particles, so is then no problem. For a coloured particle of 
width T, this is a poor approximation in the region of gluon energies (in the particle rest 
frame) below 17 The t is still sufficiently narrow that the production-decay interference 
is not a major problem, see above, but for heavier squarks and gluinos a more complex 
description may be required. The lowest coloured SUSY states, e.g. stop, tend to have 
small widths, however, and it is on such states we will concentrate our studies. The 
additional complications for very wide particles will be deferred to some future study. 

The classification by colour and spin is fairly obvious, but that does not completely 
specify the structure of the process. Consider e.g. e + e~ — > 7*/Z° — > qq. The cross section 
(neglecting mass effects) is then 

a oc e 2 c e\ + 2 e e v e KeS e q v q + {v 2 c + o^) |S| 2 (v\ + a 2 q ) , (42) 

where 

1 s 
16 sin 2 9 W cos 2 9 W s — m| + irnzTz 

represents the ratio between the Z° and the 7 propagators and couplings, 
proportional to a 2 corresponds to the qq pair coming from an axial vector 



(43) 

The term 
source, the 
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rest to it coming from a vector source. Since the QCD radiation from these two is 
somewhat different, we need to include the proper mixture, which depends on the CM 
energy, s = Eq M . It has then been simplest to perform the calculations for a pure vector 
source and a pure axial vector source, represented by 1 and 75 in Table p], and mix in the 
proportions required. We note that this mixing strategy is possible since the total cross 
section does not contain any interference terms of the character v q a q . (Such terms do 
arise when the forward-backward charge asymmetry is considered.) The pure left-handed 
mixture V — A = 1 — 75 is of special interest, since it represents the W ± bosons. Therefore 
it has been calculated separately, and is denoted by 1 ± 75 in Table [l]; as we already noted 
the sign (of the interference term in the squared matrix element) is irrelevant for the QCD 
emission aspects. In total, four alternatives are therefore open in our implementation: to 
have a pure vector source, a pure axial vector, an arbitrary mixture aV + (1 — a)A and 
the special equal mixture. In this order, that gives the four codes 11 - 14 in Table |], 
according to the numbering scheme used in the PYTHIA function PYMAEL introduced in 
version 6.153. (This routine takes as input the process code, x\, x 2 , r 2 and a, and 
returns the ratio (l/cr ) da /dxida^, omitting a factor of (a s /27r)CV.) 

Correspondingly, also most other processes can come either with or without 75 factors 
in the amplitude, or in arbitrary mixtures thereof. In the Higgs sector, normally the h° 
and H° are scalar, the A pseudoscalar and the H 1 * 1 a parameter-dependent mixture of 1 
and 75. If the coupling structure is generalized, also the neutral Higgses could be mixtures, 
however. The Qj, and qjj squark partners of the left- and right-handed quarks come with 
wave function factors 1 — 75 and 1 + 75, respectively. The squark mass eigenstates will 
be mixtures of these, with significant mixing expected especially in the third generation. 
With two squarks in a process, 7! = 1 ensures that the matrix element still can be 
written as a sum a 1 + (1 — 0)75. Again, therefore, we have chosen to perform most of 
the calculations with and without a 75 factor and then leave open to have the mixing 
depend on the current parameter choice. As a further simplification, the 1 ± 75 mixture is 
used whenever the correct choice is not known, since this mixture represents an average 
behaviour. In some instances, further restrictions exist, e.g. a pseudoscalar cannot decay 
to two scalars, at least among the MSSM processes at our disposal. 

The process at the top of Table |l] is the spin- independent eikonal answer of eq. ([18]), 
extended from the soft region where it is intended to be valid: 

1 da 

CTq dxi dx2 



The first numerator, 2{x\ + x 2 — 1 — rf — r^), here is based on an evaluation of the 2pip 2 
numerator of eq. (|1^) with p\ and p 2 given by their values after the emission of ^3. Away 
from the soft-gluon limit, however, there is some leeway in this assumption. If instead 
the pi and p 2 values before the emission of p^ had been used, one would have obtained 
x\+x 2 — l = 1— £3 — > 1. Such a seemingly minor substitution has quite dramatic effects for 
collinear emission even at rather small 23, however, and is not really an option. In order to 
have a not too unrealistic alternative to compare with, we therefore only allow deviations 
proportional to x\. (This is also what comes out of the x\ + x\ numerator of the massless 
process V — > qqg.) The extreme in this direction would be x\ + x 2 — 1 — > 1 — £3 + £3, 
but we will also allow arbitrary admixtures x\ + x 2 — 1 — > 1 — 23 + ax\, with a as a free 
parameter, obviously in no physics relation to the a introduced above. 



2{x\ + x 2 — 1 — r\ — r\ 



— C 

2n } (1 + r\ — r\ — X\){\ + r\ — r\ — x 2 ) 



2r\ 



2r\ 



'1 + r\ 



£2 



1 + rf-ri- xx. 



(44) 
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In all the formulae, the decay product mass ratios are kept as free parameters, r± = 
^bl m a and r 2 = m c /m a , while the gluon is massless. Since we are not interested in cross 
sections per se, but in the probability for gluon emission, in the LO cross sections only 
the mass dependence is retained, normalized to unity for r\ = r 2 = 0, e.g. 

a (V, A - qi q 2 ) = ~ {2 - (r? - r\f -r\- r\ ± Q ri r 2 } J (1 - r\ - rlf - *r\r\ . (45) 

What is omitted is then some set of couplings and propagators. Exactly the same set 
is also omitted from the first-order do jdx\dx 2 , thereby leaving the ratio unchanged. 
Additionally the common factor (a s /27r)Cj?, is omitted, to leave the choice of a s (p 2 ± ) free 
to be made elsewhere. The ratio of the first to leading order cross sections then gives the 
assumed differential gluon-emission rate. In the case of a mixture without and with 75 
factors, the sum has to be taken for numerator and denominator separately, e.g. 

1 da , , _ „ /rr(1 _ x ada v /dxi dx 2 + (1 - a) da A /dx 1 dx 2 

(e + e -> 7 /Z u -»• qq) = — * . (46) 



a dxi dx 2 ~ ~ a ctq + (1 — a) Oq 

The vector fraction a can here be read off from the lowest-order expression in eq. (f42|), 
where the mass factors of eq. ( |45"D were omitted. The mass effects are instead included 
in the individual da/dx\dx 2 and <r terms of eq. fl4"6|) . This standard e.g. means that 
V — A corresponds to a = 1/2, which then because of mass effects gives a somewhat larger 
vector than axial fraction. 

Note that first-order corrections to the total cross section are not included in the Uo 
denominator. This is not a unique choice, but a rather natural one: if we consider the 
gluon emission rate as the ratio of two cross sections, then including 0(a s ) corrections to 
one but not the other is not likely to improve the overall accuracy of the calculation. In 
the soft-gluon limit, it would even break the spin independence of the radiation pattern 
(cf. the next subsection), i.e. give the wrong physics. And since a complete one- loop 
calculation of both quantities is well beyond the scope of this article, we remain with 
lowest non-trivial order for both quantities. Furthermore, if the total cross section is 
written in the form cr tot = cr (l + aa s /n), then decay rates calculated so far tend to 
give a values of order unity [|28|, |29fl , i.e. small effects. An exception would be Coulomb 
corrections in the threshold region, but there the phase space for real gluon emission is 
vanishingly small anyway, so of no physical interest. By contrast, one- loop corrections to 
three-jet rates tend to be larger |BD[, although most of that is absorbed by the parton- 



shower choice of a smaller kinematics-dependent scale like p ± [30, 31 . 



4.2 Radiation patterns 

Before applying the matrix elements to specific physics situations, it is interesting to 
compare them between each other under similar conditions. This will provide an under- 
standing of the amount of spin and colour dependence in the matrix elements, and thus 
the extent to which a process-dependent Monte Carlo implementation can be expected to 
provide an improvement over a process-blind one, e.g. based on some simple dead cone 
formula. 

Fig. is intended to provide a first qualitative glimpse of differences. The amount 
of detail may be bewildering, but at this point the idea is to bring up differences and 
similarities in a broad sense, without concentrating on each process specifically. 

First consider the full curves in Figs. |^a-c, which are all for the colour structure 
1 — > 3 + 3 and the same kinematics, and only differ by the spin pattern: 1 — > 1/2 + 1/2, 
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Figure 3: The gluon emission rate as a function of the emission angle #i 3 . Specifically, the 
vertical axis gives {1 / ao)d<j / dx\ dx2 with a normalization factor (a s /27r)(5V removed. This 
three-jet phase space density differs from da/d#i3 dx% by a simple Jacobian. A variety of 
different processes are compared: decay of a colour singlet to a triplet plus an antitriplet 
full curves, ditto in the eikonal approximation dotted, decay of a triplet to a triplet plus 
singlet dashed, and gluino processes dash-dotted. The four frames differ in the scaled 
masses r\ = mi/Ecu of the two decay products, and in the gluon energy fraction x$. 
Further explanations are given in the text. 
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0^1/2 + 1/2, 1 -»• + 0, -v + 0, and 1/2 -> 1/2 + 0, in relevant cases with and 
without a 75 (but no intermediate mixtures such as 1 ± 75). In the first frame, where 
the gluon is soft, these curves almost completely overlap. At smaller x% the agreement 
becomes even better, and one can then truly speak of a universal soft-gluon emission 
pattern, with a characteristic dead cone of opening angle 9 ~ 2r = 0.4 = 23° in this 
example. By contrast, when X3 is increased, the curves tend to disagree more and more. 
At X3 = 0.3, the dead cone is still visible, although it is starting to fill in, but at 23 = 0.6 
it is completely gone. Note that the vertical scale changes significantly between the three 
frames; in absolute numbers the differences between the curves is about a factor two larger 
at £3 = 0.1 than at 0.6, but in relative terms the differences are negligible at small x%. 

This should bring home the message that the dead cone concept can only be used at 
small gluon energies, and is irrelevant for energetic gluons. It is only the lowest of the 
full curves, the completely spinless process — > + 0, that does preserve the exact dead 
cone concept, i.e. has a cross section that always vanishes in the collinear limit (in the 
rest frame of the decaying particle). 

The lower dotted curve is the eikonal expression of eq. (@), which is constructed 
to have an exact dead cone. It does agree fairly well with the spinless process, but 
not with anything else. The upper dotted curve is the modified eikonal, with an x\ 
term added. It tends to overshoot all processes at small angles, while some intermediate 
mixture, a ~ 0.5, would do a sensible job for many processes in the small-angle region. For 
medium small gluon energies and large angles, both eikonal forms undershoot, however. In 
general, one may conclude that the eikonal formula is not particularly useful for practical 
considerations, and that the process-specific matrix elements need to be used. 

Of course, the detailed pattern depends on the masses. In Fig. |3|d these are lower than 
in the first three, and representative for Z° — > bb. Then the dead cone concept works up 
to somewhat larger gluon energies, and a trace of it is still left at £3 = 0.6. It is again 
preserved exactly for — > + and the eikonal, and also approximately for 1-^0 + 0. 

In order to quantify the difference in the total amount of gluon radiation, we compare 
three measures on the three-jet phase space. Since the total three-jet cross section contains 
a soft-gluon divergence, this soft region has to be avoided. One measure is thus to integrate 
the total amount of radiated gluon energy, in shorthand 

^3 = [ x 3 — -3-^3 — dxi dx 2 ■ (47) 
J er dx 1 dx 2 

Correspondingly we define a shorthand j (1 — X\)(l — x 2 ) as an alternative removal of the 
denominator of the matrix elements (generalized to J (1 + r\ — r\ — X\){\ + rf — rf — x 2 ) 
for unequal masses). Finally a Durham distance |S[ yij = min(x^,x|)(l — cos9ij)/2 
is used to define a hard three-jet region for which = min (3/12, 2/13, 2/23) > 0.1 and a 
corresponding three-jet rate J 9(yu — 0.1), where 9 is the step function. To simplify a 
comparison between the processes, all results have arbitrarily been normalized to those 
for the V — > qq process. 

As can be seen in Table the three measures give about the same message, namely 
that differences between the processes are significant. In J X3 the ratio between the two 
extremes is a factor 1.30, in f (1 — xi)(l — x 2 ) 1.70, and in J 9(yv — 0.1) 1.91. This shows 
a steady progression of larger ratios the more one is biased towards the hard three-jet 
region. The eikonal expression is below all the calculated processes, and the modified 
eikonal (full blast) above. 

The results in Table || are for the case where the daughter masses constitute a signif- 
icant fraction of the energy available. The other limit, where instead the daughters are 
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colour 


spin 


75 


fx 3 


/(l-xi)(l-x 2 ) 


/%d-0.1) 


1^3 + 3 


1 2 ~ 2 


1 


1.000 


1.000 


1.000 






75 


1.056 


1.112 


1.133 




u 2^2 


1 


1.134 


1.293 


1.376 






75 


1.093 


1.207 


1.271 




1^0 + 


1 


1.073 


1.205 


1.310 




0^0 + 


1 


0.875 


0.758 


0.720 




1 -> i + 

2 2 ~ u 


1 


0.953 


0.918 


0.916 






75 


1.057 


1.132 


1.179 


1^3 + 3 


eikonal 


— 


0.802 


0.695 


0.659 




eikonal +£3 


— 


1.201 


1.518 


1.670 


3^3 + 1 


2 2 ~ 


1 


0.323 


0.306 


0.287 






75 


0.356 


0.365 


0.349 




1^1 + 

2 2 ~ u 


1 


0.312 


0.284 


0.258 






75 


0.357 


0.363 


0.344 




0^0 + 1 


1 


0.287 


0.242 


0.218 




0^0 + 


1 


0.279 


0.224 


0.194 




o^i + i 

w 2^2 


1 


0.359 


0.379 


0.375 






75 


0.347 


0.354 


0.346 




2 u ~ 2 


1 


0.294 


0.257 


0.239 






75 


0.314 


0.302 


0.298 


3^3 + 8 


o^i + i 

w 2^2 


1 


1.634 


1.833 


1.922 






75 


1.574 


1.712 


1.775 




i^o + i 

2 w ~ 2 


1 


1.385 


1.320 


1.291 






75 


1.549 


1.664 


1.675 


8^3 + 3 


1-^1 + 

2 2 ~ u 


1 


0.561 


0.493 


0.445 






75 


0.621 


0.607 


0.574 



Table 2: Three measures on the amount of gluon radiation in different processes, for mass 
ratios r% = = 0.2. For clarity, results have been normalized to the process in the top 
line. See the text for a detailed explanation. 



massless, is shown in Table |3[ The main message is that the process dependence remains 
also in the latter case, even if normally reduced in magnitude by about a factor of two 
between the extremes. The detailed picture is not so simple, however. Some processes 
agree in the massless limit when they do not for nonvanishing masses while, in the other 
extreme, others disagree even more for vanishing masses. 
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colour 



1^3 + 3 



3 + 3 



3 + 1 



3 + 8 



8^3 + 3 



spin 



1 

1 



1 + 1 

2 ~ 2 

+ 

+ 

1 + 



eikonal 
eikonal +x\ 



l 

2 

1 

2 






1+1 

| + 
+ 1 
+ 



+ i 



1 + 



/(l-xi)(l-x a ) 



1.000 
1.167 
1.000 
0.667 
0.917 



0.667 
1.667 



0.347 
0.347 
0.222 
0.222 
0.389 
0.264 



1.701 
1.389 



0.573 



/%D-0.1) 



1.000 
1.184 
1.141 
0.773 
0.979 



0.773 
1.595 



0.282 
0.282 
0.214 
0.214 
0.328 
0.260 



1.659 
1.384 



0.487 



Table 3: Two measures on the amount of gluon radiation in different processes, for mass 
ratios r\ = r 2 = 0, cf. Table 0. For massless daughters the 75 factor makes no difference, 
and so those results are not shown separately. The J 23 measure now is collinear divergent 
and therefore not shown. 



The set of processes with colour structure 3 — > 3 + 1 can be viewed as crossed versions 
of the 1-^3 + 3 ones, but the differences in kinematics result in another overall picture. 
Examples of radiation patterns are shown by dashed curves in Fig. ^|. The small-angle 
behaviour displays the universal dead cone effect at small X3 and again diverges wildly 
at larger 23. Note that we have oriented all processes so that the radiating daughter 
colour charge is at 0°, wherefore the radiation continues to drop off at large angles. This 
is unlike the previous set of curves, which turn around at or near the mid angle to the 
other radiating daughter. (For small £3 the two daughters are almost back-to-back, i.e. 
a bisector at 90°, while an X3 of 0.6 allows an almost symmetric configuration with 120° 
between all three.) 

It is here interesting to remind of the QED answer for the process V — > f?7 with 
rrif = mq — 



1 do- a em xl + xl ( 1-xi l-x 2 , , ., 

e f ej , (48 



a dxida;2 2n (1 - xi)(l - x 2 ) \ x 3 x 3 

where ey = ef + tj has been used to eliminate the explicit appearance of terms corre- 
sponding to radiation off the V. For ey = 0, ef = — e$ this is the QED analogue of eq. (M), 
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Figure 4: Test of the additivity of gluon emission rates, da/dx\ dx 2 with a normalization 
factor <r Cf a s /2iT removed. The dotted curves show two processes with colour flow 
3 — > 3 + 1 and 3 — > 1 + 3, respectively (and spin 1/2 — > 1/2 + 1). The upper and lower full 
curves give the complete expressions for 1 — > 3 + 3 (spin 1 — > 1/2 + 1/2) and 8^3 + 3 
(spin 1/2— > 1/2 + 0) processes, respectively. The two dashed curves, almost completely 
hidden by the full ones, are the same processes according to the additive approximations 
in eq. (|50|). 



while the current case corresponds to ej 
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da 



o« 



,-y»2 I /y»2 



a dx\ dx 2 



2n CF \l- Xl )(l-x 2 ) 



0, a 

1 — X\ 

X3 



em*3f 



F ■ 



* C _Z + A Izfl. (49) 
2tt x 3 (1-x 2 ) x 3 



Here the first part of the final expression essentially is the shower answer of eq. (^). In 
the soft-gluon limit (1 — X\)/xz ~ (l + cos6'i3)/2 (1 — cos6 l 23)/2, so this extra factor 
gives a further dampening at large emission angles, above the shower ansatz (which in 
itself is dampened by one such angular factor relative to the colour singlet decay). A 
main consequence of this extra factor is that the three-jet activity in a colour 3^3 + 1 
event is less than half of that of a 1 — >• 3 + 3 one, a pattern which remains when masses 
are included, see Table [| 

Finally we come to the processes with a gluino, denoted by dash-dotted lines in Fig. |^. 
In the limit of infinitely many colours Nc, the process 3^3 + 8 may be viewed as 
one colour flowing through from the initial triplet to the gluino, and a different colour- 
anticolour pair created between the final triplet and the gluino, i.e. as the sum of 3 — » 1+3 
and 1^3 + 3. Interference terms between the two colour flows would be suppressed by a 
factor 1/Nq, just as in the radiation pattern for the colour-related process V qqg Jl9| . 



This process therefore shows the most three-jet activity of the ones studied, especially in 
the not-displayed gluino hemisphere. 

By the same token, 8^3 + 3 may be approximated by the incoherent sum of a 
3 — ^3 + 1 and a 3 — > 1 + 3 radiation pattern, which gives less radiation than the 1^3 + 3 
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processes. The QED formula, eq. (|i5D , may here offer a convenient starting point. The 
f — > 3 + 3 process corresponds to ef = — ej, so that the interference term —e^ej is positive. 
The 8-^3 + 3 process, e.g. g — > qq, can instead be emulated in QED by ef = ej, i.e. a 
doubly-charged gluino |34j. The interference term then is negative, although suppressed 
by a colour factor 1/Nq in the QCD case. 

Thus, given a radiation pattern f 1(613) for emission off parton 1 like in a colour 3 — > 
3 + 1 process, and a corresponding fziQiz) for a colour 3 — > 1 + 3 one (essentially obtainable 
by swapping the kinematics of the above process, x\ <-> X2, ri <-» one may guess at 
the radiation pattern for the 1^3 + 3 and 8^3 + 3 processes: 

A-*+3(*is) = h<M + /apis) + 2 V/i(«i3)/a(ftLs) , 

/ 8 _ 3+ 3(^13) = /l(M+/2(0 13 )^V*)a) • (50) 

This turns out to be a good approximation up to fairly large 23 values, i.e. so long as the 
spin structure is not too important, see e.g. Fig. [|. The importance of the interference 
term also gives a simple explanation for the difference between the height of the peaks in 
the full and dashed curves in Fig. |3]a-b. 



4.3 Parity dependence 

Let us further quantify differences induced by having or not a 75 factor in the matrix 
element, i.e. between vector (V) and axial vector (A) sources, between scalar (S) and 
pseudoscalar (P) ones, etc. We will use a measure which is the mean of the matrix 
element ratios over the whole phase space for gluon emission. Since the matrix elements 
with and without the 75 coupling both have the same divergence structures the ratio is 
well-behaved everywhere. All phase space points are here given equal weight, so large 
differences may not always be reflected in significant changes of the radiation pattern. To 
be specific, we study the function 

(Q)(n,r 2 ) = r 1 — / Q(x 1 ,x 2 ,r 1 ,r 2 )dx 1 dx 2 , (51) 
J dx x dx 2 J 

where 



Q(*i, x 3> n, r 2 ) = ^ L ;z Z ■ ( 52 ) 



(Tq dcr 75 /dxida;2 
o"q 5 du 1 /da;ida;2 

X2, 0, 0) is equal to unity in all of phase space and in most cases Q(xi,x 2 , r, 0) is also 
equal to unity. The exception is when the r = mass is a boson (spin or 1) and the r ^ 
one is spin ~. All the processes in Table |I] which have both non-75 and 75 couplings are 
represented in Fig|| Processes with the same spin structure have the same ratio, where 
both the spin of the initial and final states, as well as the order of the decay products, 
are of importance. We notice that the difference between e.g. vector and axial vector 
couplings can be rather large in some cases, even for small and intermediate masses. The 
most important case, however, with two fermions in the final state, show small differences 
for r < 0.5, so this aspect will only be significant in the case of top production. Physical 
consequences will be investigated in Section |5|. 
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Figure 5: (Q)(ri,r 2 ), eq. fl5ll) , with r\ fixed, r 2 fixed, or r = r\ = r 2 . The processes are 
grouped according to spin structure. 

5 Applications 

In this section, the matrix element corrected shower is used to study some processes at 
current and future colliders. It is not to be seen as a comprehensive review, but as simple 
examples intended to illustrate the main features. 

5.1 Bottom in Z° decay 

We start by examining a process where data with good statistics already exists and 
detailed comparisons are possible, namely gluon radiation off bottom quarks produced in 
the process e + e~ — > Z° — > bb at the Z° pole. Gluons are not observed directly in the final 
state, instead they materialize as jets of hadrons. We will mainly consider jets on the 
parton level, but we also study the effects of hadronization and decays. The experimental 
data with which our results are compared have already been corrected to the parton level, 
using, among others, the same model for fragmentation which we will use. 

Jets are constructed by considering all pairs of particles in an event, where the 
particles can be hadrons, partons or clusters of hadrons or partons, finding the pair with 
the smallest "distance" y^. If this number is smaller than the "jet resolution" parameter, 
y c , the pair is joined into a new cluster by summing up their four-momenta. This procedure 
is repeated until all are larger than y c . The resulting clusters are the jets in the event 
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at the resolution scale y c . Jet algorithms differ mainly by the definition of y^, and several 



jet measures have been proposed in the literature [^. 35]. Following the lead of the LEP 
analyses we want to compare with, we have settled on the Durham [32] algorithm which 
defines 

2min(E t 2 ,Ej)(l-cos%) 

Vij = • ( 53 ) 

For small angles, this is approximately the relative transverse momentum squared of the 
pair, scaled to the visible energy in the event. If jets are constructed with a large y c , only 
the most energetic partons will be resolved, and most events will become 2-jet events. 
When y c is decreased, smaller structures will start to emerge and multi-parton final states 
become more important. It is in this region that the parton shower approach to gluon 
radiation is most useful. 

For a primary quark flavour q, the n-jet rate is defined as 

Rl(y c ) = am ~* n jcts(z/c) , (54) 

^qq— >hadrons 

which for n > 2 is generally a decreasing function of y c . We want to study the difference 
between gluon radiation off b quarks and off light quarks. A suitable observable is then 
the ratio between the respective n-jet rates 

R n(Vc) = ^H, with I = u, d or s. (55) 



Experiments at LEP |£6, [37] have found both Kg and i?| to be smaller than one, 
approaching unity for large y c . This can be understood qualitatively as a consequence of 
the well-known dead cone effect || , stating that the radiation of collinear gluons off heavy 
quarks is suppressed. The production of well separated jets, however, is not significantly 
suppressed and should approach that of the light quarks for large y c . 

We want to study this effect more quantitatively by using our improved, matrix el- 
ement corrected, parton shower. The new approach, described in Section || has been 
implemented in Pythia 6.153. This is compared to Pythia 6.152 containing the older 
approach of Section |2.4| , which includes the correct massive matrix element correction in 
the first emission only. As a reference, we also include results from the algorithm im- 
plemented in Pythia 6.129, where mass-effects are incorrectly included in the matrix 
element correction of the first emission, cf. Section |2.3| . 

The jet rate reflects both the amount of energy radiated and the direction in which it 
is radiated. Fig. |6|a shows the gluon radiation pattern in the process e + e~ — > Z° — > qq at 
the Z°-pole, for the old and the new shower routines both for light (u, d and s) quarks 
and heavy (bottom) quarks. The angle, 9 g , is defined as the angle between the radiated 
gluon and the primary quark in the CM system of the primary quark pair. In principle, 
gluons are radiated by the qq dipole as a whole, but in our implementation each radiated 
gluon is assigned to an initial quark, and in the collinear limit this separation is quite 
sensible. In the case of the total energy flow, Fig. |6|b, we add the two contributions, and 
include the energy taken by quarks produced in gluon splitting, so two symmetrical peaks 
appear. Fig. ||c shows the distribution of AE = E Tlld /EcM, i-e. the total energy fraction 
radiated in the shower. Fig. |^d, finally, shows the gluon multiplicity in the shower, which 
obviously is quite dependent on the lower shower cut-off Qq. 

We note first of all that the difference between light and heavy quarks is largest in 
the region of small angles, the dead cone, and also that the radiated energy fraction is 
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Figure 6: Gluon radiation and energy flow in the process e + e~ — > Z — > qq, where q 
is either a light (u, d, s) quark or a b quark. The b-mass is here set to the default in 
Pythia, 4.8 GeV. (a) Energy weighted gluon angle distribution, (b) Energy flow, (c) 
Radiated energy fraction, AE = E^/Eqm- (d) Gluon multiplicity. 

much larger for light quarks, peaked around 0.5 if Q = 1 GeV is used to remove the 
collinear emission off light quarks. For heavy quarks, this collinear emission is regulated 
by the mass. For large angles, the light and heavy quark distributions converge as they 
should. In the new approach, the amount of gluon radiation in b events at small angles is 
somewhat enhanced relative to the old one, where the dead cone was slightly exaggerated, 
but differences are small. By the change of evolution variable, eq. (PE[), the fraction of 
events without any shower at all (AE = 0) has decreased by almost a factor of 3 from 
about 4 per mille to 1.4. As a consequence, the peak at AE = has been flattened out. 
The gluon multiplicity has also increased slightly. 

We now want to investigate how the changes in the algorithm affect and R\ . Since 
one is studying small deviations from unity, these measures are very sensitive to changes 
in the gluon radiation pattern. Also, the LEP experiments have large samples of Z° 
events and studies at the per cent level are feasible. We will mostly study the behaviour 
of the model on the parton level, comparing different alternatives. Here the parton level 
is defined by the partonic configuration at the shower cut-off scale Qo, below which no 
further emissions occur. But first we should comment on the effects of fragmentation. 
Below the Qo scale, the Lund string fragmentation model |T] describes how the partons 
transform into the primary hadrons. Subsequently these may decay further. Once a Qo 
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Figure 7: Study of fragmentation and decay effects, (a) Effects of changing the cut-off, 
Qo, in the shower, both on the parton level and at the level of primary hadrons (i.e. 
before decays), (b) Effects of hadronization, decay of rank 1 hadrons and decay of all 
hadrons. In (b) each alternative is represented by the ± la curves given by the Monte 
Carlo statistics. 

has been chosen, the parameters of the fragmentation model should be fitted to data. Here 
we will only study the variation with Qo, both at the parton level and the hadron level, 
without any such retuning. At the hadron level, we consider first the primary hadrons 
and study decay effects separately. 

Fig. |7ja shows the effect of fragmentation and the Q dependence. We calculate Rf(y c ) 
at the parton level and at the level of primary hadrons for Qo = 1 and 2 GeV, respectively. 
The Qo dependence is largest on the parton level with i?| slightly lower for the larger Qo- 
This is because the increase in Qo increases the 3-jet rate for light quark events, while 
heavy-quark events are not influenced as much: collinear emissions, which dominate at 
the end of the cascade, are suppressed anyway. Intuitively one might have expected the jet 
rate to be higher for a lower Q , since this corresponds to a larger number of partons, but 
the further partons emitted between 1 and 2 GeV cannot give rise to new jets of their own. 
They only smear out the energy of the existing jets and possibly make them fall below 
the cut-off. Furthermore, the hadronization stage tends to decrease the Qo dependence 
because the Lund model is infrared safe and collinear gluon emissions do not affect the 
string fragmentation process. A retuning of the fragmentation parameters would further 
limit the effect of small changes in Qo- 

Fig. ^b shows the effects of fragmentation and decay. The rank 1 hadrons are the 
ones that contain the primary quarks from the decay of the Z°. In the case of bottom 
production, the rank 1 hadrons will be mainly B mesons and Ab baryons. Because of 
the hard fragmentation function for heavy hadrons, the multiplicity of primary hadrons 
in b events is smaller than that in a light quark event, reflected in the lower Rf curve. 
However, once the heavy hadrons have decayed, their decay products will more than 
compensate this. Thus, if the primary hadrons are allowed to decay, the value of R$ lies 
above unity for y c > 0.03. So, on the one hand, heavy hadrons take a large part of the 
primary quark energy, but they also decay to many particles. For there is not a large 
difference between allowing all primary hadrons to decay or only the rank 1 ones, which 
can be seen in the small difference between the two top curves in Fig. [7|b. In the following, 
we will compare our results to experimental data on the parton level, bearing in mind 
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Figure 8: R^ l (y c ) for different model parameters and variations, (a) Different main ver- 
sions. Some variations of the latest version (6.153): (b) Different bottom masses, (c) 
Different sources, (d) Other minor variations of the main theme, see the text for details. 
The data are from ]36| and [|37|. The data points from Aleph (except the first) have 
been shifted .001 units in y c for clarity. In (a) and (b) the ±lcr curves are given, while 
only the central value is shown in (c) and (d) for clarity, but the statistics is comparable. 



that this comparison can be ambiguous in view of the dependence on Qq. Our aim here 
is not to achieve a perfect fit to data, but merely ensure that the results are at the right 
level. As we will see, the data is anyway not good enough to make precise discriminations 
between all different model variations. 

In Figs. -Rg and are shown as functions of y c for several model variations. 
Figs. §a, |9|a show the difference between the old and the modified shower models. An 
even older version, where mass effects are exaggerated, is also shown as reference. Each 
model curve is displayed as a one sigma band, showing the size of the Monte Carlo error 
for 15 • 10 6 events of each kind. From the increase in multiplicity and energy flow in b 
events in the new model, one would naively expect the 3- and 4-jet rates to increase for 
heavy quarks, thus increasing R™. Actually, the main effect is to reduce the 3-jet rate 
for heavy quarks, again illustrating that allowing more collinear/soft emissions will not 
necessarily increase the rate of well separated jets. The naively expected increase of R h £ 
(from subsequent branchings) is balanced by a corresponding "smearing" loss as for R^ 1 , 
giving only a small net effect. 

As already noted, the old and new models differ mainly in the treatment of subsequent 
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Figure 9: R%(y c ) for different model parameters and variations. See caption Fig. 
details. 
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branchings; both include matrix element corrections to the first gluon emission on each 
side of the event. A clear discrimination between these two versions is not possible, 
especially in view of the variations that follow below and the relatively large experimental 
errors. It is clear, however, that the algorithm in Pythia 5.129 is ruled out by the data. 

Next we consider variations to the new default, in order to assess the uncertainties of 
the model: 

• The bottom mass. 

• Parity dependence (vector vs axial vector source). 

• Change of p\ argument in a s from z(l — z)m 2 a to z(l — z)m\(l — ml/ml) 2 in a 
branching a — > be, with a being a heavy quark. 

• Coherence effects. 

The main free parameter is the bottom mass, and Figs. |8|b, []b show the result of varying 
this mass between 4.6 GeV and 5 GeV. Fits of NLO QCD calculations to LEP data give 
a value around 3 GeV for the running bottom mass in the MS renormalization scheme 
|5|, [36|, [37] at renormalization scale mz. Our model is based on LO matrix elements and 
the mass in our case is the constituent quark mass, so they need not agree. However, a 
somewhat lower mass than the default seems to be favoured, especially for . 



In sections |4.2j and |Q| we saw that there are slight differences between different sources. 
The relevant mixture for Z° is given by eq. (|42|), but in Figs. ||c, |9]c we study the two 
extreme cases of a pure vector source and a pure axial vector one. At the Z°-pole, the 
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axial vector component dominates, as can be seen in the figure, but at higher energies the 
vector one will take over. We see that the differences are non-negligible, but not as large 
as the shown mass dependence, at least for small y c . 

Finally, in Figs. |d, |]d, some further aspects of the model, discussed in Section |3]4], 



are varied. Both the modified p± in a s and the "minimal" change of coherence tends 
to reduce R$, again showing how an increase in the total amount of radiation need not 
give more separate jets. The "intermediate" coherence option, which still is a realistic 
alternative, does give an increase, especially of R®. Thereby the overall agreement with 
data is improved, as much as with a reduced b mass, although uncertainties are sufficiently 
large that no firm conclusions should be drawn. 

Other minor issues could be studied and here we just mention a few. The definition of 
R% used here is based on a classification of events with heavy or light primary quarks. An 
alternative is to use a ratio between b-tagged and anti b-tagged events. The difference here 
is in the classification of events with gluon splitting, g — > bb, in the shower. For the R% 
ratio this is a minor issue, but it is non-negligible in the case of R?l ||. This uncertainty 
is already included in the error-bars for the experimental points in Fig. so need not be 
studied separately here. Another issue is the energy dependence of the results. Obviously 
the mass effects will decrease for larger energies, because of the decrease of r = rrib/y/s, 
and this has already been studied at 189 GeV ||. At the same time, however, the fraction 
of vector coupling increases significantly above the Z°-pole, thus decreasing R^, cf. Fig. |8|c. 
At intermediate energies this could therefore give rise to a partial cancellation of effects. 

5.2 Bottom in Higgs decay 

We now turn to gluon radiation in Higgs decay. In the standard model, the decay H° — > bb 
is expected to be large, or even dominate, for Higgs masses up to the W + W~ threshold, 
and in extensions to the standard model also heavier Higgs states can have significant bb 
branching ratios. When the fragmentation function for B mesons, i.e. the distribution 
of z = 2E-q/^/~s, is measured at LEP1, the b quark is produced from a spin 1 source. If 
instead the source is a Higgs boson, the difference in gluon radiation could give rise to a 
changed fragmentation function. Such a change would influence the experimental vertex 
detection efficiency which, if uncorrected, would give rise to errors in the determination of 
cross sections. It is therefore important to be able to describe in detail the gluon radiation 
pattern in this decay. A different measure of the nature of the source is the jet topology, 
which is also studied in the form of a modified -R3 ratio, cf. the preceding section. 

We choose to study the production of bb pairs at a CM energy of 130 GeV. This number 
is in the middle between the current lower limit and the W + W~ threshold, at around the 
value expected for the MSSM h°, but obviously the relative comparison of the sources is 
only mildly energy-dependent. Five different sources are compared. The first sample is 
considered as a reference and consists of gluon radiation in the decay 7*/Z* — > qq, where 
q is a light u, d or s quark. For simplicity they are assumed to be produced in equal 
amounts. The other four samples consist of gluon radiation in the decay X — > bb, where 
X is a vector, axial vector, scalar or pseudoscalar source. Clearly the 7*/Z* is a mixture, 
eq. (H2|), and also h°/H°/A° need not be pure states, but the separation allows us to study 



the extreme range of possibilities. 

Fig. |T0| a shows the angular energy flow. The difference between having or not having 
a 75 in the coupling is negligible, but the radiation at large angles is larger for the spin 
source than a spin 1 one. The measure in Fig. |10|b, i.e. the normalized 3-jet rate, is most 
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Figure 10: Comparison between vector, axial vector, scalar and pseudoscalar sources of bb 
pairs at 130 GeV. (a) Energy flow in the shower, (b) 3-jet rate, normalized to 7*/Z* — > qq, 
i.e. R h i{y c ) = R\{X -> hh)/R\{Y/1* -> qq), where X is V, A, S or ? and q is a light 
flavour, ler Monte Carlo error bands are shown, (c) Distribution of = 2E^/ 'y/s at the 
parton level, (d) Distribution of z B = 2E B /y/s, i.e. the fragmentation function, at the 
hadron level (only primary hadrons are considered). 

sensitive and shows that the 3-jet rate is significantly larger for a spin source, especially 
for well separated jets (large y c ), which is consistent with the larger energy flow at large 
angles. The parity dependence is very small, on the other hand, and is neglected in the 
following. 

Considering the large spin dependence one could expect the fragmentation function 
also to change considerably when going from a vector to a scalar source. Fortunately, 
as can be seen in Figs. |T0"|c-d, the changes are minor, indicating that the fragmentation 
function is mainly sensitive to the bulk of gluons at smaller angles, where the sources 
give the same emission rate. The larger gluon energy flow for a scalar source is reflected 
in a slightly smaller (z), with a difference of about 1% both at the parton and hadron 
level. The small effect on the fragmentation function is positive from an experimentalist's 
point of view, in that fragmentation functions measured at the Z° pole can be simply 
extrapolated also for possible spin sources at higher energies. 
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Figure 11: Gluon multiplicity and energy flow in the tt shower at 500 GeV centre-of-mass. 
5.3 Top production and decay 

For top production around the tt threshold, real gluon emission will be limited because of 
phase space, so aspects such as the form of the QCD confinement potential, best tackled 



with other perturbative techniques, will be more relevant p8fl . In this section we study 
the process e + e~ — > tt at 500 GeV, where gluon emission will start to become important, 
and compare the new shower routine with the older versions. We focus on fundamental 
aspects such as gluon radiation patterns and top mass determinations. In the process we 
will also study gluon radiation in the decay of the top, and differences in the radiation 
patterns for the two decay channels t — > bW + and t — ► bH + to demonstrate the advertised 
process dependence of the new shower routine. 

We first study gluon radiation in the production of the top. The matrix element 
used in the correction of the shower is an energy dependent mixture of vector and axial 
vector ones. High above the Z° pole, the vector part will dominate and at 500 GeV the a 
parameter, introduced in Section |4.1| to parameterize the relative mixture, is about 0.78 
for tt production. The dead cone effect is expected to be very pronounced in this case, 
because the scaled mass is large, around r = 0.35. Fig. |TT| shows the gluon multiplicity and 
the energy flow in the tt shower. Most events have no gluon emission, but the no-gluon rate 
has decreased in the new shower as expected. The result of an even older version, where 
the dead cone effect is severely exaggerated in all branchings, is shown as a reference. The 
lesson is that the correct description of the first emission, i.e. the difference between 6.129 
and 6.152, is much more important than that of subsequent emissions, between 6.152 and 
6.153. Once a gluon is emitted, it tends to be quite energetic and splits into several 
further gluons. This is the reason for the large spread in the multiplicity and the dip in 
this distribution for the old shower routine. The total energy flow has also increased in 
the new shower relative to the older ones. Even if the energy flow seems to be large at 
large angles, most energy is kept by the original top quarks at 6e = and n respectively. 

The dominant top decay is t — > bW + , where further gluons will be radiated. In the old 
shower routine, gluon radiation off the b did not take into account the full t-b interference 
structure. Furthermore, gluon radiation in a hypothetical decay t — > bH + was the same 
as in the one above, apart from the difference induced by the different W boson and Higgs 
masses. In order to separate off this trivial mass effect we put me = ^Wi well aware that 
this is not realistic. 
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Figure 12: (a) The gluon emission rate as a function of the emission angle #13 for top 
decay, t — > bX. The emission rate is given for the two gluon energies 23 = 0.1 (full 
curves) and 23 = 0.4 (dashed). The curves corresponding to the higher energy have been 
multiplied by a factor of 40 for clarity. Three different curves are shown for each energy: 
pseudoscalar Higgs, W boson and scalar Higgs (from high to low curves), (b) The 3-jet 
rate in complete tt events as a function of y c for the three different cases above, where 
the W/H are removed from the cluster search, (c-d) Gluon multiplicity and energy flow 
in the shower induced by the top decay (in the rest-frame of the decaying top with the 
W at 9 E = tt). 

The matrix element used in the decay t — ► bW + has the familiar V — A coupling, 
but the t — ► bH + case is less obvious. In the MSSM the vertex factor is proportional 
to m t (l — 75) + mt,tan 2 /?(l + 75), i.e. a parameter-dependent mixture of scalar (S) 
and pseudoscalar (P) couplings. Naively, the S — P one would be expected to dominate 
because of the large difference between the top and bottom masses. However, the favoured 
values of tan (3 > 3 could cancel this effect and give rise to an almost pure scalar coupling. 
We therefore study the two extreme cases of S and P separately. The difference between 
the matrix elements of this type has been studied in Tables ^| and |3] and found to be 
non- negligible for massive daughters but to vanish in the massless limit. For the process 
t — > bX, where X is a W + , pure scalar H + or pure pseudoscalar H + , some differences are 
visible in the gluon energy flow for larger gluon energies, see Fig. [L^a, but the differences 
induced in observables such as the jet rates are small, Fig. |T2]b. That is, we are closer 
to the massless case of Table [| indicating that the small mass of the radiating b is more 
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Figure 13: tt events after top decays, (a) Total gluon multiplicity, (b) The 3-jet rate on 
the parton level, with W/H removed from the clustering. Initial-state photon radiation 
is included, but both models are equally affected by this, (c) Reconstructed top mass, 
parton level, (d) Reconstructed top mass, hadron level. 



important than the large mass of the non-radiating W + /H + . In what follows, therefore, 
only the decay to W + is shown. This is the most important case and, as it turns out, the 
intermediate alternative, cf. Fig. pk. 



Figs. |l2|c-d show the gluon multiplicity and energy flow in the decay of the top for 
the new and the older shower routines. The gluon multiplicity is slightly more peaked 
in the new shower and the no-gluon rate has been somewhat reduced. The angle 9e is 
defined as the angle of emitted energy to the primary bottom quark in the rest frame 
of the decaying top, which is why no peak appears at large angles, corresponding to the 
direction of the colour neutral W/H. In the new correction to the shower, gluon radiation 
at large angles is more severely suppressed. As a consequence, there is more energy left 
for radiation at small ones, allowing the curves to cross there. This is an effect in addition 
to the influence of the new matrix element correction scheme, visible e.g. in Fig. |6|a. The 
difference between decays to W and H are very small also here (not shown). The result 
of Pythia 6.129 is similar to 6.152, but the dead cone is slightly more pronounced, as 
noted before. 

In a complete event, gluon emission from the decay and production are added, and 
also possible radiation in hadronic decays of the W/H must be considered. In principle, 
these showers will also interfere but, as a realistic first approximation, we work in the zero 
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width limit where these interferences vanish. Considering the two stages of production 
and decay together, but ignoring the decay products of the W/H, Fig. |T3| a shows the total 
gluon multiplicity for the new and older shower routines. Since the radiation has increased 
in the production and decreased in the decay, effects tend to cancel. The net result is a 
slight increase in the average gluon multiplicity and a somewhat smaller width. The 3-jet 
rate, however, has decreased, Fig. |13|b, presumably because of the reduced probability to 



have gluon emissions at large angles in the top decay, cf. Fig. |12|d. 

The total gluon multiplicity is not a very good measure of the event properties and 
certainly not an observable. To give a practical example where gluon radiation is impor- 
tant, we consider the reconstruction of the top mass in a simplified scenario, assuming 
that the W/H can be completely reconstructed in order only to study the effects of gluon 
emission in the top production and decay. We find two jets and calculate the top mass by 
considering the two possible combinations of jet+W momenta. The combination which 
minimizes (mi — m t ) 2 + [m-i — m t ) 2 , where m t = 175 GeV, gives the reconstructed top 
masses mi and m^. The distribution of the reconstructed mass is given in Fig. p~3] c, where 
jets are clustered on the parton level. The reduced radiation level survives in this measure 
as a slightly narrower peak for the new shower. The effects of fragmentation and decay 
are quite large for the reconstructed mass, but some differences still survive in the high 
mass wing, as seen in Fig. |T3|d. 



5.4 Supersymmetry production and decay 

In supersymmetric models, a whole new set of processes involving coloured particles is in- 
troduced. We have calculated the most important LO matrix elements of the MSSM, listed 
in Table [I], to be used as input to the parton shower routine. In versions of Pythia prior 
to 6.153, supersymmetric particles did not shower. Because of the large masses of these 
hypothetical particles, this is a good first approximation. However, if supersymmetric 
particles are discovered at the LHC or a future linear collider, detailed studies of their 
properties will profit from a better understanding of QCD radiation patterns. 

A complication that appears in the MSSM is significant rates of three body decays 
involving one or several coloured particles in the initial or final states. Often several 
interfering diagrams with intermediate off-shell propagators contribute. An example is 
the decay x + ~ > x'Vq with either a W or squark propagator. These processes do not 
fit into the present framework with corrections to matrix elements of the type a — > beg. 
As a preliminary solution, the q'q pair above is assumed to be produced from a V — A 
source. A new kind of gluon radiation process in the MSSM is from four-vertices such as 
q — > q'W + g. These diagrams introduce no new divergences. The interference terms only 
contribute to the ordinary divergences, and the squared amplitude gives rise to a constant. 
An upper estimate is again given by the new shower expression, which is corrected by the 
full matrix element. 

As a first example, consider top vs. stop production at a linear e + e~ collider and the 
decay t — ► bW + vs. t — ► bx + . Again, identical masses are assumed and consequently 
mf = m t and m^ = m\y are used. As usual, the tt source is an energy dependent 
mixture of vector and axial vector. Only a scalar coupling is possible for Z° — > tt, 
so no ambiguity here either. The top decay has the V — A coupling, while the stop 
one is a parameter dependent mixture of scalar and pseudoscalar. Again we consider 
the two extreme cases separately. While the stop obviously decays isotropically, the 
current Pythia implementation does not include top polarization information and thus 
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Figure 14: Top vs. stop production and decay at a 500 GeV e + e collider, (a) energy 
flow in the tt (tt) shower, (b) 3-jet rate on the parton level for the production and decay 
taken together, with the W/x removed from the clustering. 



also decays this particle isotropically. Initial-state photon radiation has not been included 
in the comparison, since the different threshold behaviours of the tt/ tt cross sections allow 
more energetic radiation in the former process, which will affect event topologies. 



Fig. T^a shows the energy flow for top vs. stop production. The multiplicity is the 
same for both top and stop, but the energy flow is slightly larger for stop at large angles. 
The difference between the gluon radiation patterns in the top and stop decays is negligible 
(not shown). As a result of the increased gluonic energy flow, the 3-jet rate at the parton 
level, with the W/x removed from the cluster analysis, is slightly larger for stop, Fig. [TJ|b. 
We notice further that the parity dependence is negligible also in this case. The differences 
between supersymmetric and standard model processes is not large in this example, at 
least when mass effects are neglected, and the ambiguity in the coupling structure is 
negligible. This is all good news, but since gluon emission from supersymmetric particles 
is a completely new feature, many more tests should be done. 

As a last example we consider a full simulation of gluino pair production (via gg — > 
gg) at a 14 TeV pp collider. Since we have not developed a full shower formalism for 
2^2 processes, as a first approximation the eikonal expression will be used for radiation 
off the gg system, including a colour factor rescaling by Nq/Cf = 9/4. Initial state 
radiation is studied separately and multiple interactions are neglected. A scenario with 
m g = 450 GeV, = 250 GeV and = 200 GeV is used. The dominant gluino decays 
are then g — > bib (tit) and the squarks decay by processes of the type q — > xb- If the x is 
a neutralino, which here is the lightest supersymmetric particle, it will not decay further. 
A chargino decays to a neutralino plus either leptons or quarks in approximately equal 
amounts. This MSSM scenario is implemented in the Pythia event generator |J and the 
details are described in the Spythia manual |39| . 



To assess the importance of gluon radiation off supersymmetric particles in a hadron- 
hadron collision, we compare the transverse gluon energy flow from SUSY and non-SUSY 
particles in the final state radiation (FSR), Fig. |l5|a. The additional gluon radiation 
from coloured SUSY particles is thus seen to be small compared to the ordinary gluon 
radiation from quarks. In this paper we have only considered QCD final state radiation 
in e + e~ annihilation. In a hadron collider environment, gluons can be radiated also from 
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Figure 15: Gluino production at a 14 TeV pp collider, (a-b) Transverse gluon energy 
flow in pseudorapidity, r] = |ln((p +p z )/(p — Pz))- ( c ~d) Jet multiplicity for R = 0.75, 
E± t min = 10 and 40 GeV respectively. 

the incoming partons (initial state radiation, ISR). This contribution is almost as large 
as the FSR one, Fig. |15|a, but the rapidity distribution is different. 

In Fig. [I5| b the three Pythia versions considered in this paper are compared. Only 
the FSR component is shown, and the additional radiation from SUSY particles in the 
latest version gives rise to an increased total transverse gluon energy flow. This increase is, 
however, slightly compensated by the decrease in gluon radiation in the top and squark 
decays, cf. Section |5.3| . In Figs. |I"5|c-d the result of a simple cluster search is shown. 
Particles with a summed transverse energy, ^2 E± > E^ min inside a a/ (Arj) 2 + (A</?) 2 < R 
cone are joined in a cluster. For large E^ min , the jet multiplicity is slightly increased 
when gluon radiation from supersymmetric particles is included. This could be caused 
by energetic gluons radiated off the gluinos. For small Ej_ tmin , on the other hand, smaller 
structures are probed and the radiation off b and lighter quarks is most important. In 
summary, we conclude that gluon radiation off supersymmetric particles at a hadron 
collider is small compared to the other sources, and of importance mainly in high-precision 
studies. 
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coherence 


g -> uu + dd + ss (%) 


g - cc (%) 


g - bb (%) 


<*b(P±) 


full 


14.3 


1.26 


0.15 


<*b(P±) 


intermediate 


14.8 


1.27 


0.16 


<Xb(P±) 


reduced 


21.1 


1.92 


0.26 




none 


38.8 


3.06 


0.31 


a s (m 2 /4) 


full 


12.9 


1.15 


0.15 


a s (m 2 g /4) 


intermediate 


13.3 


1.17 


0.15 


a s (m 2 g /A) 


reduced 


20.0 


1.78 


0.28 


a s (m 2 /4) 


none 


43.3 


3.47 


0.46 



Table 4: The rate of gluon splitting to qq pairs, in Z° decay at 91.2 GeV with the normal 
primary flavour mixture. 



6 Gluon splitting to heavy quarks 

Data at LEP1 show a larger rate of secondary charm and bottom production than pre- 
dicted in most shower descriptions || |40[1 , or in analytical studies [41|. We therefore 
comment on a few of the issues in the Monte Carlo simulation and how a relaxation of 
some demands would affect rates. 



6.1 Strong coupling argument and kinematics 

The default behaviour in Pythia is to let a s have p\ as argument. Actually, the ex- 
act kinematics has not yet been reconstructed when a s is invoked, so the approximate 
expression p 2 ± z(l — z)m^ is used, see discussion at eq. (|7|). Since a s blows up when 
its argument approaches Aqcd, this translates into a requirement on p 2 ± or on z and m a , 
restricting allowed emissions to p± > Qo/2, where Q is the shower cut-off scale. Also 
when full kinematics is reconstructed, this is reflected in a suppression of branchings with 
small p±. Therefore, in g — > qq branchings analyzed in the g rest frame, the quarks do not 
come out with the 1 + cos 2 9 angular distribution (with respect to the direction of motion 
of the gluon) one might expect away from threshold, or a somewhat more isotropic one 
closer to threshold, but are rather peaked at 90° and dying out at 0° and 180°. 

For g —>■ qq branchings, the soft-gluon results that lead to the choice of p\ as scale 
42| are no longer compelling, however. One could instead use some other scale that does 



not depend on z but only on m a = m g , i.e the off-shellness of the branching gluon, and 
remove the pj_ cut. A reasonable choice, even if not unique, is to use m 2 /4, where the 
factor 1/4 ensures continuity with p\ for z = 1/2. This possibility has been added as a 
new option. 

Actually, the change of a s argument in itself leads to a reduced g — > qq splitting rate, 
while the removal of the p± > Qo/2 requirement increases it. The net result is a decrease 
by about 10%, Table [|. The topologies of the events are changed somewhat, so rates 
within experimental cuts could be more affected. However, the changes are not as big as 
might have been expected — see the following. 
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6.2 Coherence 



In the above subsection, it appears as if the 1 + cos 2 9 distribution would be recovered 
in the new a s (m 2 /4) option. However, this neglects the coherence condition, which is 
imposed as a requirement in the shower that successive opening angles in branchings 
become smaller. Such a condition actually disfavours branchings with z close to or 
1, since the opening angle becomes large in this limit, eq. (|8]). It should be noted that 
the opening angle discussed here is not the true one, but the one based on approximate 
kinematics, including neglect of masses. More generally, the coherence formalism is not 
really developed with this kind of configurations in mind, especially not with a heavy 
quark pair close to threshold. 

As a means to exploring consequences, two new coherence level options have thus 
been introduced. In the first, the p\ of a g — » qq branching is reduced by the correct 
mass-dependent factor, 1 — 4m 2 /m 2 , while the massless approximation is kept for the 
longitudinal momentum. This is fully within the uncertainty of the game, and no less 
reasonable than the default. In the second, no angular ordering at all is imposed on 
g — > qq branchings. This is certainly an extreme scenario, and should be viewed with 
caution. However, it is still interesting to see what it leads to. 

In Section [3]4] another "intermediate" coherence variation of the default "full" one 
was introduced, affecting the rate of gluon emission off the primary quarks but not the 
subsequent gluon cascades. This variation has negligible consequences for the secondary 
quark rate at LEP energies and is not considered further. 

It turns out that the decay angle distribution of the gluon is much more distorted 
by the coherence than by the a s and kinematics considerations described earlier. Both 
modifications are required if one would like to have a 1 + cos 2 9 shape, however. Also 
other distributions, like gluon mass and energy, are affected by the choice of options. 

The most dramatic effect appears in the total gluon branching rate, however, Table |4j. 
Already the reduced angular ordering requirement can boost the g — > bb rate by almost 
a factor of two. The effects are even bigger without any angular ordering constraints at 
all. It is difficult to know what to make of these big effects. Experimental information on 
the angular distribution of secondary cc/bb pairs might help understand what is going 
on, but probably that is not possible experimentally. Anyway, the measured (but still 
uncertain) values are fully bracketed by the range of the models, indicating that there 
need not be a conflict between theory and experiment. 



7 Summary and Outlook 

We have in this article studied QCD radiation off heavy particles, both based on the 
calculation of a wide set of first-order matrix elements and by the usage of these matrix 
elements as input to a parton-shower description of multiple-gluon emission. 

The matrix-element calculations provide at least two important insights. One is the 
significant spin dependence on the rate of gluon emission, given identical kinematical 
conditions. Some such effects could be expected simply from the different shape of the 
splitting kernel of a quark and a squark, say, but equally important is the spin of the 
decaying particle, and the parity of the process (combining that of mother and daughters). 
The best illustration is the difference between a vector and a scalar colour singlet source 
each decaying to a pair of squarks, where the gluon emission rate differs by up to a factor 
of 1.8 in our limited study, and is likely to reach even higher levels in other observables. 
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The 75 dependence disappears in the massless limit, but the bulk of other effects remain 
also there, so this is not only an issue for the production of heavy flavours. This process 
dependence would not be caught by the traditional parton-shower philosophy, where the 
universal aspects of gluon emission are emphasized. 

The other insight is that the 'dead cone' concept is one that only applies universally 
for soft gluon emission. Only one process provides an exact zero in the collinear emission 
rate while, for reasonably hard gluons, some processes show no dip at all at small angles. 
And, while soft gluons may be ideologically interesting, the experimentally observable 
event shapes are often more crucially dependent on some intermediate range of gluon 
energies. Again, therefore, the message is that a universal shower description may be 
misleading. 

Fortunately, it is not impossible to combine the shower picture with a process depen- 
dence of the kind noted above. We have in this article developed one approach to the 
problem, wherein matrix element information is provided for all shower branchings of the 
primary particles. This combines a process dependent resummation of gluon emission off 
these particles with the traditional strengths of the shower formalism, such as providing 
exclusive final states with exact energy-momentum conservation. 

To illustrate this formalism, a few physics examples have been studied. For bottom 
production at the Z° pole, detailed comparisons with data are possible. We find the new 
shower routine to be in good agreement with current data if the uncertainties in both 
model and data are taken into account. A significant dependence on coherence options 
were found both for jet rates and gluon splitting rates in the shower. It has further been 
possible to rule out an older version of the shower routine where mass effects were not 
correctly accounted for. 

Differences in jet rates were found for different coupling structures in Higgs decay to 
bottom, whereas the influence on the fragmentation function is minor. Both the produc- 
tion and decay of the top quark is dependent on the matrix element correction, influencing 
e.g. the mass reconstruction of the top quark. Again the jet rate was found to be most 
sensitive. 

Gluon radiation off supersymmetric particles, squarks and gluinos, has been intro- 
duced, but the effect of this additional radiation is small in high energy processes, es- 
pecially at hadron colliders, because of the background from initial state radiation and 
showers off the standard model decay products. Differences between top and stop events 
are small if equal masses are assumed, but many more comparisons between standard 
model and supersymmetric processes could be envisaged. 

From this limited study we conclude that gluon emission off b and lighter quarks 
dominate the picture. If this component is modelled well, the rest is less important, but 
still of ideological interest. Since the scaled bottom mass r = m h /m x , in the decay of a 
resonance X, is small in most cases, we see little if any parity dependence. In Z° decays 
this effect is less than 1%, but it could be visible in a high-precision study. Also for the 
spin dependence, effects are small in most cases, although the H° — > bb example shows 
they need not be quite as small. It thus turns out that the new process dependence will 
mainly be important when high precision is strived at. 

We note that the current study is not the end of the story. While we have sampled a 
fair selection of colour and spin structures, more are likely to turn up even in the simple 
context of two-body decays considered here. The MSSM also allows a significant rate 
for three-body decays, normally as a sequence of two consecutive two-body ones with an 
intermediate off-shell propagator, and with the possibility of interference between several 
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such intermediate states. The real challenge, however, may well be provided by the more 
complex production processes at hadron colliders, where the concept of a sequence of 
s-channel processes need no longer be valid, e.g. squark production gg — ► qq with a t- 
channel squark propagator. The large rate of initial-state radiation also implies that the 
initial-final state interference terms may be more important than the final-state radiation 
off the heavy particles themselves. Much work therefore lies ahead, if one desires a good 
description of QCD effects in many processes of interest at the LHC. 
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